AFRL-PR-WP-TR-1999-2078 


STUDIES  OF  SURFACE  DEACTIVATION 
OF  VIBRATIONALLY-EXCITED 
HOMONUCLEAR  MOLECULES  IN 
GASEOUS  DISCHARGE  MEDIA  USING 
COHERENT  ANTI-STOKES  RAMAN 
SPECTROSCOPY  (CARS) 


PERRY  P.  YANEY,  PH.D. 
JOHN  W.  PARISH 


UNIVERSITY  OF  DAYTON 
300  COLLEGE  PARK 
DAYTON,  OH  45469 


JANUARY  1999 


FINAL  REPORT  FOR  04/01/1993  -  08/01/1998 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


199910U  OU 

PROPULSION  DIRECTORATE 
AIR  FORCE  RESEARCH  LABORATORY 
AIR  FORCE  MATERIEL  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE  OH  45433-7251 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose  other  than 
in  connection  with  a  definitely  Government-related  procurement,  the  United  States  Government 
incurs  no  responsibility  or  any  obligation  whatsoever.  The  fact  that  the  government  may  have 
formulated  or  in  any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be 
regarded  by  implication,  or  otherwise  in  any  manner  constmed,  as  licensing  the  holder,  or  any 
other  person  or  corporation;  or  as  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell 
any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  is  releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it 
will  be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


Power  Systems  Branch 
Power  Division 
Propulsion  Directorate 


HORWATH 


Acting  Chief,  Power  Systems  Branch 
Power  Division 
Propulsion  Directorate 


iR^YEBEAM 
'Acting  Chief 
Power  Division 
Propulsion  Directorate 


If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list,  or  if  the 
addressee  is  no  longer  employed  by  your  organization  please  notify  AFRL/PRPS, 
Wright-Patterson  AFB,  OH  45433-7919  to  help  us  maintain  a  current  mailing  list. 

Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security  considerations,  . 
contractual  obligations,  or  notice  on  a  specific  document. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704^0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  Instructions,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing 
the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information 
Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  {0704-0188),  Washington,  DC  20503. 


1.  AGENCY  USE  ONLY  {Leave  blank}  2.  REPORT  DATE 

JANUARY  1999 


4.  TITLE  AND  SUBTITLE 

STUDIES  OF  SURFACE  DEACTIVATION  OF  VIBRATIONALLY-EXCITED 
HOMONUCLEAR  MOLECULES  IN  GASEOUS  DISCHARGE  MEDIA  USING 
COHERENT  ANTI-STOKES  RAMAN  SPECTROSCOPY  (CARS) _ 


6.  AUTHOR(S) 

PERRY  P.  YANEY,  PH.D. 
JOHN  W.  PARISH 


3.  REPORT  TYPE  AND  DATES  COVERED 


FINAL  REPORT  FOR  04/01/1993  -  08/01/1998 


5.  FUNDING  NUMBERS 
C  F33615-93-C-2303 
PE  61102 
PR  2301 
TA  PI 
WU  01 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

UNIVERSITY  OF  DAYTON 
300  COLLEGE  PARK 
DAYTON,  OH  45469 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

UDROOOlZ 


9.  SPONSORINGIMONITORING  AGENCY  NAMEIS)  AND  ADDRESS(ES) 

PROPULSION  DIRECTORATE 

AIR  FORCE  RESEARCH  LABORATORY 

AIR  FORCE  MATERIEL  COMMAND 

WRIGHT-PATTERSON  AFB,  OH  45433-725 1 

POC:  RTSWA  N.  GANGin.Y.  AFRL/PRPS.  937-255-2923 


12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 


APPROVED  FOR  PUBLIC  RELEASE,  DISTRIBUTION  UNLIMITED. 


13.  ABSTRACT  (Maximum 200  mrds} 


SEE  REVERSE 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

AFRL-PR-WP-TR-1999-2078 


14.  SUBJECT  TERMS 


Nitrogen,  Wall  relaxation  rate,  Heterogeneous  deactivation  coefficient,  Coherent  anti-Stokes 
raman  spectroscopy,  CARS,  Vibrational  deactivation  rates  of  nitrogen  X  state 


115.  NUMBER  OF  PAGES 


116.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


20.  LIMITATION  OF 
ABSTRACT 


Standard  Form  298  (Rev.  2-89)  (EG) 

Presented  by  ANSI  Std.  239.18 

Designed  using  Perform  Pro,  WHS/DIOR,  Oct  94 


ABSTRACT 


Deactivation  coefficients,  Yv»  of  vibrationally-excited  NjCX^Sg^v)  on  alloys  of  aluminum, 
stainless  steel  alloys,  a  titanium  alloy,  gold,  Pyrex  glass  and  Teflon  were  determined  from 
measurements  of  the  wall  deactivation  rate,  k,,  for  the  v  =  1  to  4  or  5  hot  bands  using 
coherent  anti-Stokes  Raman  spectroscopy  (CARS).  A  water-cooled  positive  column 
discharge  was  used  to  populate  the  vibrational  states.  The  nitrogen  gas  flow  rate  was 
adjusted  to  optimize  the  observed  decay  of  the  populations.  The  excited  gas  was  exposed  to 
the  tubular  specimen  surface  for  precisely  controlled  residence  times  and  efforts  were  made 
to  track  the  influence  of  surface  properties.  The  BOX  CARS  configuration  was  used  to 
spatially  resolve  only  the  gas  that  was  exposed  to  the  specimen  surface.  Extensive 
measurements  on  Pyrex  gave  Yi  values  between  2.4x10'^  and  6.7x10'^  depending  on  the 
treatment  history  of  the  surface.  The  values  for  Y4  ranged  from  2. 9x  10"^  for  the  AMS  4943D 
alloy  of  titanium  to  approximately  unity  for  the  AMS  312  stainless  steel  alloy.  The  low  value 
for  titanium  can  be  attributed  to  the  oxide  layer.  The  observed  values  for  gold  decreased  with 
increasing  temperature  consistent  with  a  physical  adsorption  process.  For  nearly  all 
specimens,  k^  was  proportional  to  v,  or  nearly  so.  Moreover,  Yv  values  for  neat  Pyrex,  gold 
foil  and  aluminum  showed  similar  proportional  dependencies  on  v. 
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INTRODUCTION 


Vibrationally  excited  molecules  play  various  roles  in  gas  kinetics  and  plasma 
chemistry.  Of  the  seven  diatomic,  homonuclear  molecules  occurring  at  normal  temperatures, 
nitrogen  holds  interest  due  to  its  importance  in  combustion,  biology,  microelectronics,  high 
speed  and  high  altitude  aircraft  and  spacecraft  reentry.  Vibrational  energy  in  the  ground 
electronic  state  of  nitrogen  is  important  in  the  short-lived  pink  afterglow  production  where 
Penning  ionization  and  the  rapid  e-V  exchange  leads  to  considerable  heat  loss  and  N2(B^Zu^) 
pumping.'  Vibrational  energy  is  important  in  promoting  plasma  chemistry  where  it  is  an 
available  energy  source  for  molecular  dissociation.  Here  radicals  or  atoms  are  delivered 
where  needed  as  in  plasma-enhanced  chemical  vapor  deposition  (PECVD).^  Published 
research  into  the  catalytic  chemistry  of  active  nitrogen  in  low  pressure  post-discharge  gas 
revealed  the  presence  of  a  large  concentration  of  excited  ground  state  nitrogen  NjCX'Sg^,  v) 
accompanying  N(‘*S3/2)  atoms.^  Accurate  calorimetric  measurements  of  atom  recombination 
require  correction  for  the  vibrational  energy  deposition.  Accurate  theoretical  models  of 
vibrational  distribution  require  specific  state-selective  measurements,  which  have  been 
developed."*’^ 

Vibrational  energy  transfer  ofN2(X'Sg‘^,  v)  to  solid  surfaces  is  addressed  in  this  work. 
Population  distributions  and  loss  rates  in  the  presence  of  a  surface  were  studied.  These 
results  extend  the  knowledge  about  diatomic  molecular  wall  loss  rates  to  higher  v-states  than 
is  currently  available.  This  information  can  influence  the  design  of  gas  lasers  and  improve  our 
understanding  of  all  reactive  gaseous  environments,  such  as  in  combustors,  in  plasma 
processing  chambers  for  microelectronics  and  on  high  temperature  exterior  surfaces  of  high 
speed  aircraft.  In  this  work,  vibrationally  excited  nitrogen  was  exposed  to  tubular  reactor 
surfaces  of  common  construction  materials.  This  list  includes  neat  Pyrex,  Pyrex  exposed  to 
active  nitrogen,  type  3003  aluminum,  t)q>es  304  and  AMS  312  stainless  steel,  AMS  4943D 
titanium,  99.999%  pure  gold  and  Teflon.  The  tubular  reactor  elements  had  an  internal 
diameter  of  about  1.5  cm  and  were  cut  into  specimens  about  12  cm  long. 
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A  flowing  afterglow  created  in  a  water-cooled,  dc  electric  discharge  reacted  with  each 
surface  for  precisely-stepped  residence  times.  The  glow  discharge  provided  a  low 
translational  temperature  and  a  high  vibrational  excitation  with  insignificant  atom  and  charged 
species  densities.  This  gave  spectra  which  had  isolated  vibrational  bands  with  significant 
rotational  line  intensities  limited  to  approximately  J  =  21 .  The  vibrational  wall  deactivation 
rates  were  measured  for  the  v  =  1  to  either  the  v  =  4  or  5  vibrational  “hot”  bands  by  probing 
the  axial  flow  with  a  small-volume  optical  probe  using  coherent  anti-Stokes  Raman 
spectroscopy  (CARS).  Although  signals  from  the  v  =  6  band  were  observed,  useful 
measurements  above  the  v  =  5  band  were  virtually  impossible  under  our  conditions.  In 
essence  this  means  those  populations  were  too  small  to  have  an  impact  on  the  experiment  as 
a  whole.  Nonequilibrium  distributions  were  characterized  by  rotational  temperatures  of  ~3  50 
K  and  vibrational  temperatures  of  -2000  K.  Numerous  experiments  were  made  on  some 
surfaces  to  determine  repeatability.  At  the  same  time,  a  slow  change  in  the  surface  due  to  the 
afterglow  was  observed  for  every  surface.  Measurements  were  made  on  the  axis  of  the  gas 
flow  and  referred  to  the  surface  by  calculation. 

Observed  deactivation  rate  k„  which  include  both  homogeneous  and  heterogeneous 
deactivation,  and  the  “apparent”  deactivation  coefficients  for  the  surfaces  are  reported  for 
each  measured  v  v  +  1  band.  The  vibrational  distribution  in  the  bulk  of  the  gas  depends 
on  the  deactivation  rates  and  radial  diffusion.  Relaxed  molecules  must  diffuse  to  the  probe 
volume  before  the  effect  can  be  seen.  Moreover,  the  observed  deactivation  due  to  the  surface 
is  significantly  modified  by  the  fast  vibrational-vibrational  (V-V)  exchange  processes  via 
diffusion.  These  processes  soften  the  effect  of  wall  deactivation  by  the  redistribution  of  the 
populations  up  the  vibrational  energy  level  ladder.  This  process  is  strong  when  there  is 
significant  populations  in  the  low-lying  vibrational  levels  such  as  produced  by  a  gaseous 
discharge.  For  the  low  lying  v-levels  studied  here  (^5),  the  vibrational-translational  (V-T) 
exchange  rates  are  insignificant.  In  order  to  extract  the  “true”  wall  deactivation  rates  k^^) 
from  k„  a  Master  Equation  Model  of  Nj  was  fitted  to  the  data  for  four  specimens  with  the 
k^v^)  as  fitting  parameters.  This  model  includes  V-V  and  V-T  energy  exchange  processes 
as  well  as  wall  exchange  interactions.  It  uses  the  Schwartz- Slawsky-Herzfeld  (SSH)  law  of 
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scaling  deactivation  rates  and  the  development  of  the  model  closely  follows  that  of  Bailey® 
and  Capitelli/  From  this  procedure,  “true”  values  of  and  were  determined.  As 
described  in  a  later  section,  this  fitting  effort  was  greatly  complicated  by  the  fact  that  the 
model  computed  populations  out  to  v  =  40  while  data  was  available  only  up  to  v  =  4  or  5. 
This  required  an  extrapolation  procedure  to  constrain  the  fitting  of  the  populations  above  v 
=  4  or  5.  It  was  found  that  many  of  the  data  sets  could  not  be  adequately  constrained  to 
permit  unambiguous  fittings.  Nevertheless,  it  is  suggested  that  the  four  successful  fittings 
provide  information  that  can  be  used  to  provide  reasonable  corrections  of  the  for  all  the 

specimens  to  viable  estimates  of  the  “true”  y^ 

The  Experimental  Details  section  describes  the  gaseous  discharge  system  and  the 
procedures  used  to  establish  a  clean  gaseous  system.  The  CARS  optical  system  and  the 
electronic  data  acquisition  system  are  also  described.  The  Theory  section  reviews  the  basic 
characteristics  of  the  CARS  technique  and  gives  a  detailed  analysis  of  the  dependence  of  the 
deactivation  coefBcient  measurements  on  the  observed  deactivation  rate  and  diffusion.  A 
review  of  the  Master  Equation  Model  is  also  included.  The  data  collection  procedures,  the 
corrections  for  stimulated  Raman  scattering  and  the  normalization  procedures  are  presented 
in  the  Analysis  section.  This  section  also  includes  the  procedure  used  to  measure  the  peak 
values  of  the  rotational  lines  and  the  method  use  to  unfold  the  population  decay  rates  of  the 
individual  v  states  from  the  population  differences  obtained  form  the  CARS  data.  The  Data 
and  Results  section  reports  the  results  of  all  the  deactivation  measurements  including  the 
fitting  results  using  the  Master  Equation  Model.  The  results  of  the  XPS  measurements  on 
selected  specimens  are  reported  here  also. 
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EXPERIMENTAL  DETAILS 


In  an  experimental  apparatus  similar  to  that  used  by  Black  et  al.,*  the  vibrationally- 
excited  molecules,  generated  in  the  low  pressure  dc  discharge,  were  carried  in  the  afterglow 
through  a  1 5-cm-long  tubular  reactor  in  which  diffusion  carried  them  to  the  surface.  Black 
et  al.  used  mostly  a  thermal  source  and  a  1 00-cm-long  reactor.  Because  of  this  long  flow  tube 
length,  they  had  to  correct  for  Poiseuille  flow.  Moreover,  they  could  only  achieve  sufiBcient 
excitation  to  observe  the  v  =  1  0  transition  using  spontaneous  anti-Stokes  Raman 

scatttering.  A  fi'action  y,,  of  the  excited  molecules  which  collide  with  a  surface  having  a 
catalytic  capability  result  in  the  transfer  of  vibrational  energy  to  the  surface  vibrational 
modes.®  The  population  distribution  remaining  in  the  gas  was  measured  after  stepped 
residence  times  of  contact  of  the  gas  with  the  surface.  The  pressure  chosen  for  these 
measurements  was  typically  between  10  and  20  Torr.  This  pressure  was  sufficient  to  provide 
adequate  CARS  signals  and  low  enough  to  permit  efficient  steady  state  production  of 
vibrationally-excited  nitrogen  in  the  discharge.  From  the  observed  decay  of  the  vibrational 
populations  with  time,  the  total  rate  coefficients  k^  were  determined  from  which  the  apparent 
wall  deactivation  efficiencies  were  calculated  for  each  vibrational  level.  The  total  rate 
coefficients  are  independent  of  the  size  of  the  probe  volume  as  well  as  the  radial  location  in 
the  tubular  specimen.  They  do  depend,  however,  on  the  combined  effects  of  collisions  with 
the  surface  and  the  gas-gas  collisions.  For  the  conditions  of  these  studies,  the  latter  is 
dominated  by  V-V  processes.  For  given  values  of  “true”  yy,  the  observed  rate  coefficients 

ky  depend  on  the  temperature,  pressure  and  tube  radius.  The  yy,  however,  are  characteristics 
of  the  surface  which  are  determined  by  surface  activity,  temperature,  the  colliding  molecule, 
and  possibly  pressure.  By  fitting  the  Master  Equation  Model  to  these  population  data, 
extraction  of  k^yan)  of  homogeneous  effects  was  attempted.  However,  because  the 
coupled  equations  are  stiff  and  the  fitting  process  is  highly  dependent  on  the  fact  that  the  rate 
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of  change  of  population  in  one  level  is  highly  coupled  to  neighboring  populations  via  the  V-V 
processes,  this  effort  gave  viable  results  only  for  certain  data  sets. 

A  flowing  gas  handling  system  carried  nitrogen  (Air  Products  ultra-pure  carrier  grade 
99.999%  pure)  through  a  water-cooled  dc  electric  discharge  operating  in  the  range  of  1 1  to 
22  Td  with  40  cm  of  the  positive  column  exposed  to  the  flow.  The  flowing  afterglow  passed 
around  a  right-angle  bend  and  then  through  the  tubular  reactor  specimen.  Tylan  General  FC- 
260  flow  controllers  maintained  steady  mass  flow  to  in  the  range  of  0  to  2000  seem  calibrated 
to  within  0.2%.  The  gas  supply  system  was  all  stainless  steel  except  for  the  Pyrex  discharge 
tube,  the  source  tube  (which  carried  the  active  nitrogen  to  the  specimen  tube)  and  a  0.8  m 
section  of  flexible  polyethylene  tubing  at  the  input  of  the  discharge  tube  which  prevented  the 
plasma  from  striking  up-stream.  The  experiments  were  accomplished  between  1 1  and  23  Torr 
(67  to  3066  Pa)  with  most  measurements  at  14  Torr.  In  Fig.  1,  the  positive  column  discharge 
pumps  the  vibrational  states  to  nonequilibrium.  The  vibrational  temperature  of  the  gas 
entering  the  specimen  tube,  T„  ranged  from  about  1 800  to  2800  K  while  the  gas  temperature, 
Tg,  was  between  about  3 1 5  to  380  K.  Either  a  dc  discharge  or  a  cw  microwave  source  could 
have  been  used  to  preferentially  pump  N2(X,v).“’  The  dc  discharge  avoided  high 
concentrations  of  N(‘‘S)  atoms  and  N2  B-  and  A-state  molecules.  Recombination  of  the 
residual  atom  concentration  produced  the  weak  straw-colored  Lewis-Rayleigh  afterglow.*  " 
A  spectroscopic  study  of  this  N2(B  -  A)  glow  with  residence  time  provided  a  measure  of 
atom  recombination  on  the  wall  and  in  the  gas  and  gave  a  rough  estimate  of  N  density.  *^  The 
mole  fraction  of  N  atoms  was  found  to  be  about  0.0055.  The  transit  time  of  the  gas  in  the 
source  tube  was  no  less  than  20  ms.  This  insured  that  the  density  of  electronically-excited 
molecules  was  at  least  two  orders  of  magnitude  less  than  the  N  atom  density  and  thus  a  small 
fraction  of  the  observable  N2(X,v)  populations  of  the  gas  impacting  the  specimen  surface. 
Moreover,  all  charged  species  had  sufficient  time  to  become  neutralized  well  before  they  came 
in  contact  with  the  specimen  surface. 

Exiting  the  discharge,  the  gas  turned  90  degrees,  past  a  Wood's  horn,  into  the  source 
tube  where  a  pink  afterglow  terminated.  The  line-of-sight  from  the  discharge  tube  terminated 
into  the  Wood’s  horn  in  order  to  prevent  near-UV  emissions  from  the  plasma  exciting  the  gas 
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Figure  1 .  Diagram  of  discharge  tube  and  probe  chamber  for  measuring  deactivation  rates 
as  a  function  of  Z. 

I 


entering  the  reactor  specimen  tube.  During  the  transit  to  the  specimen  tube,  the  gas  cooled 
and  attained  steady  state  conditions  giving  a  steady  population  at  the  exit  of  the  source  tube. 
The  position  of  the  source  tube,  which  was  sealed  through  an  0-ring  gland  in  the  top  of  the 
optically-accessible  main  chamber,  could  be  positioned  to  ±0. 1  mm,  providing  very  accurate 
control  of  the  residence  time  of  the  gas  in  the  specimen  tube.  Brewster-angle  windows  were 
mounted  on  15-cm-long  tubular  extensions  on  the  chamber  which  permitted  focusing  of  the 
high  peak-power  CARS  beams  in  the  chamber  without  damaging  the  windows.  There  was 
also  a  pair  of  side  windows  for  emission  experiments.  The  CARS  signals  were  generated  by 
probing  the  gas  along  a  diameter  through  the  axis  of  the  specimen  tube  in  a  1  -cm  gap  between 
bottom  of  the  specimen  tube  and  the  top  of  the  gas  “catch”  tube.  The  catch  tube  was 
connected  to  a  vacuum  pump  through  a  set  of  three  valves  arranged  to  provide  for  on/off, 
coarse  and  fine  control  of  the  pumping  rate  on  the  main  chamber.  Pressure  was  regulated  by 
manual  butterfly  and  needle  valves  some  150  cm  distance  after  the  specimen  tube  and 
monitored  with  a  pair  of  Wallace  and  Tieman  gauges  (0  to  20  and  0  to  200  Torr  absolute) 
connected  to  the  main  chamber.  Additional  6-mm-diameter  gas  outlets  in  the  floor  of  the 
main  chamber  constrained  recirculation  of  the  gas.  A  computer-generated  analog  signal 
controlled  the  mass  flow  controller.  Thermocouple  and  Convectron  gauges  were 
incorporated  for  base  pressure  and  helium  leak  checking  of  the  system,  respectively. 

At  more  than  a  meter  upstream  of  the  discharge  tube,  a  mixing  chamber  that  had  three 
inlet  ports  was  used  to  mix  nitrogen  with  argon  or  helium  as  needed.  The  system  was 
returned  to  atmospheric  pressure  with  pure  nitrogen  or  argon.  Before  each  experiment,  each 
specimen  was  degreased  with  methylene  chloride  in  a  sonic  cleaner  before  insertion  in  the 
main  chamber.  The  source  and  specimen  surfaces  were  cleaned  with  an  argon  discharge 
around  2  Torr  at  60  mA.  The  pressure  was  varied  between  0.25  and  10  Torr  during  this 
procedure  to  promote  out-gassing  of  any  trapped  gases.  No  oxygen  or  water  could  be 
detected  in  gas  samples  taken  from  the  system  and  analyzed  locally  by  Fourier  transform  mass 
spectrometry.*^  The  supply  and  exhaust  systems  were  periodically  baked  out  assisted  by  a 
combination  of  heat  tapes  and  heat  gun  with  either  flowing  gas  or  pumping  at  the  1  pm  base 
pressure  for  six  to  eight  hours.  Base  pressure  was  accomplished  with  an  isolated  secondary 
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system  including  an  Edwards  mechanical  pump  fitted  with  a  molecular  sieve  at  the  inlet  port. 
The  primary  flow  system  was  pumped  by  a  Welch  mechanical  pump  fitted  with  an  oil  mist 
filter  and  vented  outside  the  lab.  At  no  time  were  primary  and  secondary  pumps  operated 
together.  Moreover,  the  primary  pump  was  never  opened  at  or  below  its  base  pressure  of  25 
pm  and  was  always  operated  with  flowing  gas. 

Mass  flow  was  regulated  to  a  fixed  value  for  each  experiment  between  450  and  1200 
seem  ±0.2%.  This  provided  an  axial  transport  speed  between  about  2  and  8  m/s  in  the 
specimen  tube.  This  gave  from  15  to  60  ms  maximum  residence  time  of  the  vibrationally 
active  gas  in  the  specimen  tube  and  more  than  20  ms  in  the  source  tube.  Because  the 
specimen  tube  length-to-diameter  ratio  was  small  (^8),  plug  flow  was  assumed,  thereby 
requiring  no  correction  for  Poiseuille  flow. 

CARS  is  an  effective  probe  of  ground  electronic  state  vibrational  levels  in  homopolar 
diatomic  molecules.  Lack  of  a  vibrational-rotational  absorption  spectrum  due  to  the  absence 
of  electric  dipole  transitions  compels  the  use  of  a  Raman  probe.  At  low  pressures  the 
nonlinear  process  produces  little  signal  generation  advantage  over  linear  Raman  spectroscopy; 
however,  the  nonlinear  process  generates  a  signal  that  is  confined  to  a  coherent  beam  which 
greatly  simplifies  the  signal  collection.  On  the  other  hand,  two  complications  demand  careful 
attention.  The  change  in  the  third  order,  nonlinear  susceptibility  in  the  CARS  process  is 
proportional  to  the  difference  between  the  populations  of  the  two  quantum  levels  probed  by 
the  laser  beams.  This  requires  an  analysis  step,  discussed  in  the  Analysis  section,  that  extracts 
fractional  populations  from  signal  data  containing  population  differences.  Secondly, 
stimulated  Raman  scattering  (SRS)  simultaneously  pumps  population  from  the  lower  to  the 
upper  level,  which  introduces  a  laser-power-dependent  systematic  shift  in  the  signal 
measurements.  This  “saturation”  effect  has  been  the  subject  of  an  investigation  similar  to 
studies  reported  here.  However,  it  was  necessary  to  independently  determine  the  correction 
for  the  intermediate  peak  laser  power  levels  used  in  these  studies.’®  This  correction  requires 
a  calibration  step  which  is  briefly  described  in  the  Analysis  section. 

The  pulsed  lasers  used  in  the  CARS  system  consisted  of  an  injection  seeded,  doubled 
Spectra-Physics  DCR-3  Nd:  YAG  laser  which  pumped  a  Lumonics  HyperDye  300  tunable  dye 
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laser.  The  YAG  laser  produces  a  7  ns  pulse  and  provided  a  CARS  pump  beam  at  532  nm  with 
a  single  longitudinal  mode  linewidth  of  0.006  cm*’  (0.00017  nm,  180  MHZ).  The  dye  laser, 
using  sulforhodamine  640  dye,  produced  a  5  ns,  multimode  CARS  probe  beam  with  a 
linewidth  of  0.1  cm-'(0.0037nm,  3  GHz)overthe  range  of  16,640  to  16,390  cm-' (601  to610 
nm).  It  was  found  that  22.33  mg  dye  per  liter  methanol  dilution  provided  a  slowly  varying 
power  curve  over  this  range.  Figure  2  is  a  diagram  of  the  optical  system  used  in  this  work. 
The  figure  shows  how  the  appropriate  timing  delay  of  the  YAG  beams  was  accomplished 
through  a  series  of  prisms.  The  YAG  and  dye  laser  beam  energies  were  adjusted  using  half¬ 
wave  plates  in  combination  with  polarizing  prisms.'*  The  main  pump  beam  was  split  into  two 
equal-energy,  parallel  beams  about  2  cm  apart  with  a  beam-splitter  mirror  and  aligned 
symmetrically  on  the  focusing  lens.  The  probe  beam  was  increased  in  diameter  through  a 
beam-expanding  telescope  in  order  to  approximately  match  the  dye  beam  waist  to  the  waist 
sizes  of  the  two  pump  beams  in  the  focal  volume  of  the  lens.  The  probe  beam  was  aligned  in 
a  plane  perpendicular  to  the  plane  of  the  pump  beams  so  as  to  satisfy  the  BOX  CARS  phase 
matching  condition.'’  The  pump  beams  were  elliptical  in  cross  section  with  major  and  minor 
axis  dimensions  at  full  width  at  half  maximum  (FWHM)  of 300  pm  and  100  pm,  respectively. 
The  waist  diameter  of  the  dye  beam  was  70  pm  FWHM.  These  dimensions  were  determined 
by  knife-edge  test  and  bum  patterns.  The  length  of  the  overlap  or  probe  volume  of  the  three 
beams  was  measured  to  be  about  5  mm  by  translating  a  thin  argon  sheet  through  the  volume 
while  measuring  the  air  signal.  The  tolerance  on  the  positioning  of  the  probe  volume  was  ±2 
to  3  mm  due  to  the  length  of  the  probe  volume  and  the  1 5-mm  tube  diameter.  This  positioning 
was  accomplished  by  translating  the  main  chamber  parallel  to  the  probe  beams  while  measuring 
the  CARS  signal  for  the  v  =  1  to  2  transition.  Beam  alignment  was  checked  before  each 
experiment  by  peaking  the  ~473-nm  CARS  signal  with  the  dye  beam  focusing  lens.  During 
chamber  maintenance,  beam  crossing  alignment  was  achieved  by  passing  the  three  attenuated 
beams  through  a  pin  hole.  The  BOX  CARS  beam  arrangement  was  chosen  because  it 
permitted  precise  location  of  the  measurement  volume.  A  scheme  using  60“  rotated 
polarization  in  each  beam,  available  for  suppressing  the  nonresonant  background,  was  not  used 
because  the  signal  is  reduced  by  one  third.'*  Since  the  nonresonant  background  in  these 
measurements  was  very  weak,  this  signal  reduction  was  deemed  too  great  considering  the 
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Figure  2.  CARS  optical  system  used  in  measuring  the  deactivation  rates.  Labels:  AP,  aperture; 

BD,  beam  dump;  BS#,  beam  splitter;  DHS,  dichroic  harmonic  separator;  L#,  lens;  M#,  mirror; 
P#,  prism;  PBS#,  polarizing  beam  splitter;  PIN,  photodiode;  PMT,  photomultiplier;  SHG,  second 
harmonic  generator;  WP#,  half-wave  plate. 


diminutive  size  of  the  signals  from  the  low  pressure  gas.  All  beams  were  horizontally 
polarized  parallel  with  the  optical  table. 

Beam  powers  incident  on  the  main  chamber  were  tracked  with  PIN  diodes  calibrated 
between  each  run  within  each  experiment.  Diode  counts  linearly  tracked  the  incident  fluence 
and  systematic  energy  deviations  in  the  laser  pulses.  While  the  diode  signal  counts  were 
proportional  to  the  incident  pulse  energy,  average  power  measurements  were  more 
convenient  to  use  because  the  average  diode  signal  was  used  in  normalizing  the  average 
CARS  signal.  These  diodes  were  an  important  nexus  between  saturation  calibration  and  the 
deactivation  experiments.  The  count-to-mW  ratios  of  each  deactivation  experiment  and  the 
associated  saturation  calibration  experiment  were  used  to  convert  the  power  in  the 
deactivation  measurements  to  the  equivalent  power  used  in  the  saturation  measurements.  The 
degree  of  saturation  S  was  then  calculated  and  returned  to  correct  the  raw  CARS  signal 
values. 

Instead  of  using  a  high  pressure  argon  gas  cell  to  reference  the  measurements,  which 
was  really  impractical  for  our  system  besides  being  notoriously  difficult  to  keep  aligned,  the 
averages  of  the  outputs  of  the  laser  PIN  diodes  were  used  to  normalize  the  CARS  signals. 
These  integrated  signals  tracked  power  changes  in  the  two  beams  within  a  run  as  well  as 
between  runs  within  an  experiment.  During  any  given  deactivation  or  saturation  run,  the 
standard  deviation  of  the  mean  of  the  output  of  the  YAG  laser  over  200  pulses  was  well 
under  1%.  Due  to  mode  beating  in  the  dye  laser,  there  was  no  significant  correlation  between 
the  fluctuations  in  the  dye  laser  output  and  the  YAG  laser  output.  Thus,  assuming  that  there 
was  no  significant  correlation  between  the  susceptibility  fluctuations  and  the  dye  laser 
fluctuations,  the  normalization  of  the  CARS  signal  was  accomplished  by  dividing  the  average 
CARS  signal  by  the  individual  averages  of  the  laser  diode  signals.  The  normalization  and 
saturation  correction  are  detailed  in  the  Analysis  section.  In  addition,  signal  changes  due  to 
the  dye  spectral  power  curve  were  removed  initially  by  calibrating  the  wavelength  response 
of  the  dye  laser  PIN  diode.  It  was  found  subsequently  that  the  combined  spectral  response 
of  the  receiver  optics,  spectrometer  and  photomultiplier  tube  (PMT)  was  flat  over  the 
operating  range  and  needed  no  correction. 


11 


The  receiver  optics  consisted  primarily  of  a  500-mm-focal-length  planoconvex  lens, 
a  right-angle  prism  for  spatial  separation  of  signal  and  laser  beams,  two  p-polarization 
dichroic  mirrors,  two  p-polarization  dichroic  filters,  and  a  39-mm,  best-form  lens  for  imaging 
the  CARS  beam  on  the  entrance  slit  of  a  1-m  Jarrell- Ash  grating  spectrometer.  The 
spectrometer  was  used  primarily  as  an  additional  filter  for  rejection  of  any  pump  beam  leakage 
light.  An  EMI  9816KB  PMT  performed  CARS  signal  detection  typically  operating  at  2100 
V.  The  rejected  beams  were  directed  to  beam  dumps  except  for  about  4%  reflected  from  the 
hypotenuse  of  two  prisms  toward  the  PIN  diodes. 

Data  acquisition  and  dye  laser  tuning  was  computer  controlled  through  a  Keithley 
Data  Acquisition  System  Model  570.  The  data  acquisition  and  control  (DAqCt)  software  is 
written  in  GW  BASIC  which  incorporates  Keithley  SoftSOO  BASIC  libraries  for  controlling 
the  Keithley  570.  The  final  DAqCt  program  version,  named  SB9.BAS,  consists  of  three 
separate  programs:  SB9,  SB9KL1  and  SB9KL2.  This  software  has  several  features  including 
the  capability  for  manually  controlling  the  dye  laser  and  the  flow  controllers,  automatically 
tuning  the  dye  laser  from  a  script  file  and  carrying  out  the  acquisition  of  a  continuous 
spectrum  or  a  skip  spectrum  from  a  script  file. 

Gated  detection  of  the  PMT  and  PIN  diode  signal  pulses  was  used  to  minimize 
background  signals  and  noise  from  ambient  light  and  from  radio  frequency  interference 
generated  primarily  in  the  YAG  flash  lamp  power  supply.  The  system  trigger  pulse  was 
generated  by  a  PIN  diode  from  the  YAG  pulses  which  leaked  through  a  turning  mirror 
following  the  beam  splitter  for  the  dye  laser.  This  trigger  pulse  was  sent  to  a  LeCroy  Model 
821  Discriminator  to  generate  a  NIM  standard  pulse  which  triggered  a  Hewlett-Packard 
(H-P)  Model  3  0 1 2B  pulse  generator.  The  output  from  this  generator  was  attenuated  and  split 
into  a  clock  start  pulse  and  three  -1.6-V  gate  pulses,  each  70  ns  long.  The  gate  pulses 
provided  the  gating  of  the  LeCroy  Model  2249SG  fast  pulse  analog-to-digital  converter 
(ADC).  The  ADC  digitized  the  integrated  signal  pulses  from  the  PMT  (i.e.,  the  CARS  signal) 
and  the  PIN  diodes  for  the  YAG  and  dye  lasers.  These  signal  pulses  were  appropriately 
delayed  with  coax  delay  cables  such  that  they  arrived  at  the  ADC  within  the  gate  pulses. 
Each  pulse  was  passed  first  through  a  10-step  Phillips  Scientific  Model  771  pulse  amplifier. 
The  gain  was  set  to  maximize  the  utilization  of  the  ADC’s  1 000-count  dynamic  range  so  that 
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the  maximum  pulse  did  not  exceed  the  dynamic  range.  The  CARS  signal  was  additionally 
attenuated  using  a  motor-controlled  neutral  density  filter  wheel  mounted  in  fi'ont  of  the  PMT. 
This  was  necessary  because  of  the  more  than  1000  to  1  range  of  the  CARS  signals.  Both  the 
amplifiers  and  the  neutral  density  filters  were  calibrated  to  accuracies  better  than  0. 1%.  The 
YAG  laser  pulse  noise  level  of  1%  relative  standard  deviation  was  the  lowest  of  the  three 
signals.  The  dye  laser  presented  about  4%  to  6%  noise  while  the  noise  on  the  CARS  signal 
was  typically  between  22  to  33%  relative  standard  deviation.  Each  laser  PIN  diode  had  an 
adjustable  optical  attenuator  that  greatly  alleviated  spatial  alignment  requirements  and 
removed  channel  spectra  in  the  dye  diode.  The  attenuators  consisted  of  a  pair  of  telescoping 
brass  tubes  with  an  opal  glass  window  at  the  interior  end  of  each  tube. 

A  fifth  pulse  from  the  first  H-P  pulse  generator  triggered  a  second  H-P  3012B  pulse 
generator  which  produced  a  2.3  ms  pulse  for  gating  the  Keithley  570.  The  computer  reads 
the  ADC  after  each  laser  pulse  as  many  times  as  was  programmed  for  that  deactivation  data 
point.  After  each  data  point,  the  program  totaled  the  data  and  stored  the  value  in  a  file.  The 
wavelength  was  stepped  and  the  process  repeated.  Later,  the  average  signal  values  and 
uncertainties  were  then  calculated. 

The  water  cooling  of  the  nitrogen  glow  discharge  tube  kept  the  gas  temperature  at  the 
specimen  tube  below  400  K.  This  discharge  tube  was  connected  to  a  manually-operated,  high 
voltage  dc  power  supply  through  a  high-power  50  kQ  ballast  resistor.  Typical  voltages  were 
about  8.8  kV  at  the  power  supply  giving  about  4.3  kV  at  the  discharge  terminals.  The 
discharge  current  was  set  to  90  mA  and  monitored  periodically  throughout  the  experiment. 
Hollow  Kovar  electrodes  were  used  in  the  discharge  tube. 
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THEORY 


CARS  Probe  Technique 

The  CARS  technique  was  originally  developed  for  the  combustion  community. 
Combustion  is  often  accompanied  by  fluorescence  emissions,  high  luminous  emissions  due 
to  soot  formation,  high  temperatures,  chemical  reactivity,  turbulence  and  sometimes  high 
transport  velocities.  Physical  probes  can  disturb  the  chemistry  or  the  population  distributions, 
they  are  not  generally  state  selective  and  they  may  be  unable  to  physically  survive  the 
environment.  An  optical  probe  can  overcome  many  of  these  difficulties.  CARS  has  the 
highly-desired  property  of  producing  a  large,  coherent  optical  signal  beam  on  the  short 
wavelength  side  (i.e.,  blue  side)  of  the  incident  laser  wavelength  thereby  removing  the  signal 
from  the  strong,  potentially  interfering  fluorescence  and  incandescence  emissions  that  are  on 
the  red  side  of  the  laser  wavelength.  These  emissions  are  typically  generated  when  a  laser 
beam  is  incident  on  a  complex  combustion  medium.  This  is  in  contrast  to  the  spontaneous 
Raman  or  laser-induced  fluorescence  techniques  wherein  the  signals  are  usually  on  the  red 
side  of  the  laser  wavelength.  These  techniques  also  have  the  disadvantage  of  requiring  signal 
detection  over  a  significant  solid  angle  which  introduces  spatial  resolution  restrictions  and 
limitations  of  use  on  highly  luminous  media.  CARS  has  the  disadvantages  of  requiring  two 
optical  ports  that  can  withstand  high  peak  power,  short  optical  pulses  and  a  limiting  sensitivity 
of  roughly  one  part  per  thousand  (1  O’®  cm'^)  in  atmospheric  flames.  Under  the  conditions  of 
the  experiments  reported  here,  the  number  density  was  almost  two  orders  of  magnitude  below 
the  value  at  atmospheric  pressure  and  ranged  over  nearly  3.5  orders  of  magnitude  which 
brings  the  lower  limit  to  no  more  than  10'*  cm‘^.  There  are  several  appropriately  named 
variants  of  the  technique  depending  on  the  geometry  of  the  three  input  beams.  A  review  of 
these  is  given  in  Eckbreth.’*  General  reviews  of  the  CARS  technique  can  also  be  found  in 
Arnold  et  al.,”  Tolies,^®  Druet^’  and  Attal-Tretout.^^  The  essential  characteristics  of  the 
technique  are  presented  here.^ 

The  CARS  process  is  governed  by  a  fourth-order  susceptibility  tensor,  and  thus  is  a 
member  of  the  nonlinear  family  of  four- wave  mixing  processes.  The  polarization  P(t)  of  a 
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scattering  medium  is  related  to  an  incident  electric  field  E(t)  by  the  relation  P(t)  =  Sq  x  ^(0-> 
where  x  is  the  electric  susceptibility  of  the  medium  and  Bq  is  the  permittivity  of  free  space. 
Expanding  the  polarization  in  a  power  series  one  obtains 

Pit)  =  (1) 

where  the  x^"^  factors  are  tensors  of  order  n  +  1 .  In  four-wave  mixing,  the  total  applied  field 
is  made  up  of  three  incident  fields  that  can  be  written  as 

Eit)  =  (2) 

where  is  one  of  the  incident  plane  waves  of  frequency  o);  and  the  propagation 

directions  are  parallel.  All  three  terms  in  Eq.  (1 )  can,  in  general,  be  simultaneously  nonzero. 
In  a  gas,  however,  the  x^’^  term  describes  Rayleigh  scattering  of  the  incident  waves  and 
vanishes  identically.  Substituting  Eq.  (2)  into  the  x‘^’  term  of  Eq.  (1),  the  general  expression 
for  the  four-wave  mixing  polarization  can  be  written 

jkl  imno) 


(3) 

where  (mno)  indicates  the  sum  over  the  frequency  components  is  restricted  to  the  condition 
indicated  in  the  polarization  argument.  The  number  of  separate  elements  in  x^^^yu  is  8 1 .  The 
number  of  elements  that  are  actually  needed  is  determined  by  symmetry  properties  of  the 
probed  medium.  For  an  isotropic  rnedium  like  a  gas,  there  are  21  nonzero  elements  that 
possess  the  property  that  any  Cartesian  component  which  appears  at  least  once  must  appear 
an  even  number  of  times.  Thus,  because  of  the  equivalence  of  the  three  coordinate  axes,  the 
only  surviving  elements  are  xxxx,  xxyy,  xyxy,  xyyx  where  xxxx  =  xxyy  +  xyxy  +  xyyx,  which 
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specifies  the  polarizations  of  the  four  fields.  In  the  experiments  reported  here,  the 
polarizations  of  the  incident  fields  were  parallel  corresponding  to  the  xxxx  component. 

The  terms  of  Eq.  (3)  of  interest  here  are  of  the  type  given  in  Eq.  (4),  where  the  first 
expression  corresponds  to  the  case  where  there  are  two  different  pump  waves  at  ©i  and  ©2 
with  the  probe  wave  at  ©3.  The  second  and  third  expressions  are  essentially  the  same  with 
©2  =  ©1.  (Note  that  the  V2  factors  are  not  shown  on  the  right  sides  of  these  expressions.) 

P(0)4  =  c0i+a)2-W3)  = 

/>((04  =  2(i),-a)3)  =  (4) 

P((04  =  20)2-0)3)  =  3x^^^E^E* 

In  the  case  where  two  incident  frequencies,  one  for  a  pump,  oOp  =  ©,  or  ©2 ,  and  the  other  for 
the  probe,  ©,  =  ©3  (Stokes),  corresponds  to  the  resonance  condition  of  ©^  =  ©p  -  ©, ,  where 
©V  is  the  transition  energy  between  two  vibrational-rotational  levels  of  nitrogen,  then  the 
number  of  terms  contributing  to  the  same  interaction  reduces  to  the  12  terms  summarized  in 
Eq.  (4).  If  the  expression  for  the  electric  field  is  written  with  both  the  temporal  and  spatial 
parts,  where  ^  are  the  propagation  vectors,  then 

E((xi.,r)  =  ^Ef^e\p(-i(o.t)exp(ik.-r)  +  c.c.  (5) 

Substituting  into  Eq.  (4)  with  ©2  =  ©1  =  ©p,  the  spatial  phase  of  the  coherent  scattered  wave, 
called  the  anti-Stokes  wave,  at  frequency  ©4  =  ©a,  becomes 

k'f=(^k^-kyf  . 

as  ^  p  s' 

If  the  incident  waves  are  planar  with  and  parallel  and  directed  in  the  z  direction,  the 

signal  wave  will  experience  growth  in  this  direction  at  the  frequency  ©a^  =  2©p  -  ©,. 
Substituting  Eqs.  (4)  and  (5)  into  the  wave  equation,  applying  the  slowly  varying  wave 
approximation  and  defining  the  phase  mismatch  Ak  =  2kp  -  k,  -  ka^^  one  obtains 


16 


(6) 


2ik. 


dEi(Oa,,z)  w; 


c" 


as  ^2'^  -  '■  P'  ' 


Integrating  Eq.  (6)  over  z  =  0  to  squaring  the  result  to  obtain  intensity  and  using 


/,  =  ^|£(o>,)|^ 


(7) 


for  each  field,  where  n^  are  the  refractive  indices  and  c  is  the  speed  of  light,  the  expression  for 
the  CARS  signal  beam  intensity  can  be  written 
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where  Xcars  is  the  third  order  susceptibility  for  the  CARS  process.  Maximum  signal  is 
achieved  when  the  phase  matching  condition  Ak  =  0  is  satisfied.  It  should  be  pointed  out  that 
the  pump  intensity  is  on  the  order  of  1  GW/cm^  in  practice  and  the  probe  beam  is  weaker  by 
a  factor  of  roughly  10  but  provides  the  coupling  to  the  resonance  properties  of  the  medium. 
The  strength  of  the  CARS  signal  for  given  peak  intensities  of  the  pump  and  probe  beams 
depends  on  the  square  of  the  CARS  susceptibility,  Xcars  >  which  is  proportional  to  the 
spontaneous  Raman  cross  section  times  the  difference  in  the  population  densities  of  the  two 
levels  connected  by  the  process. 

Diffusion  and  the  Wall  Loss  Flux 

The  probability  of  losing  vibration  quanta  from  nitrogen  molecules  to  a  surface  was 
determined  by  measuring  the  rate  of  decay  of  the  vibrational  state  populations  in  the  center 
of  the  reactor  tube.  At  a  pressure  of  1 7  Torr  and  a  temperature  of 390  K,  the  V-V  exchange 
rate  is  about  2x  lOVs  which  is  roughly  one  tenth  the  rate  molecules  in  the  probe  volume  visit 
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the  surface  and  about  10^  times  the  wall  deactivation  rate.  There  is,  therefore,  significant 
probability  that  molecules  deactivated  at  the  wall  will  be  re-excited  by  the  V-V  process  before 
returning  to  the  axis  where  the  measurements  are  made.  Thus,  except  for  very  small 
deactivation  probabilities,  the  observed  rates  of  decay  of  the  populations  versus  residence 
time  (i.e,  the  length  of  the  specimen  tube  exposed  to  the  gas  divided  by  the  flow  velocity)  will 
be  slower  than  the  actual  wall  rates.  There  are  several  likely  mechanisms  for  direct  energy 
exchange  with  the  surface,  each  possibly  at  a  different  rate.  Some  mechanisms  result  fi-om 
an  event  sequence  referred  to  as  a  channel  wherein  there  may  be  interactions  between 
channels.  There  have  been  numerous  discussions  of  the  chemistry  and  thermodynamics  on 
surfaces  .  The  measurements  reported  here  were  incapable  of  distinguishing 

among  these  mechanisms,  but  rather  they  provided  an  integrated  evaluation  of  all  active 
surface  loss  and  gain  mechanisms.  Nevertheless,  the  accuracy  of  the  measurements  of  these 
net  loss  rates  was  not  diminished  by  this  fact. 

Molecules  are  reflected  from  the  wall  without  loss  of  vibrational  energy  with  a 

probability  of  1  -  Jy.  With  diffusion,  as  Jy  increases,  a  radial  density  profile  in  the  vibrational 
level  populations  develops.  In  the  experiments,  the  length-to-diameter  ratio  of  the  reactor 
tube  was  about  8,  which,  with  the  modest  gas  flow  rates  used,  maintained  a  flat,  “plug-type” 
radial  velocity  profile  along  the  axis  of  the  tube.  The  relationship  between  wall  loss  rate  ^v<wa//) 
and  probability  Jy  is  influenced  by  the  cylindrical  geometry  of  the  reactor  tube.  Because  the 
rate  constant  could  be  quite  high,  especially  for  v  =  4  and  5,  the  following  development 
retains  all  corrections.  (Note  that  the  “v(wall)”  label  is  dropped  temporarily.)  The  total  flux 
of  particles  (particles  per  unit  area  per  unit  time)  near  a  surface,  J,  is  given  by^* 

J  =  ,9, 

where  [A^2]  is  the  number  density  of  molecules,  v„  is  the  velocity  normal  to  the  wall  (i.e, in  the 
radial  direction)  and  ffvj  is  the  probability  density  of  velocity.  Assuming  the  Maxwell- 
Boltzmann  velocity  distribution  function,  the  flux  becomes 
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where  kg  is  Boltzmann's  constant,  Tis  the  absolute  temperature,  and  c  is  the  mean  molecular 
speed.  This  expression  must  be  modified  when  there  are  reflection  and  loss  fluxes  at  the 
surface.  The  expression  now  is  given  in  two  equations  which  describe  the  transport  to  and 
from  the  wall.  Thus, 

jf  =  (II) 

for  vibrational  species  /,  where  the  flux  toward  the  wall  is  X,  J  is  the  flux  in  the  opposing 
direction,  D  is  the  coefficient  of  diffusion  and  Wg  is  mean  transport  drift  (flow)  velocity 
towards  the  surface.  Since  the  flow  is  only  along  the  tube  axis,  Wg  =  0.  Accounting  for  the 
loss  at  the  surface,  we  have  for  that  boundary  condition 

=  (12) 

Substituting  Eqs.  (1 1)  in  to  Eq.  (12)  and  simplifying  gives 

_  r[Arjc-  _ 

■  4(1-1)  <>*) 

Equation  (13)  shows  that  for  y  =  0  the  flux  reduces  to  that  given  by  Eq.  (10).  In  situations 
where  y  is  small,  the(l  -y/2)'’  inEq.  (13)  can  be  neglected,  but  this  is  not  generally  applicable 
for  the  measurements  reported  here.  The  steady  state  diffusion  equation  can  be  written^'* 
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dz 


where  Mq  is  the  axial  transport  velocity.  A  solution  of  the  form  [Ay  (r,  z)  =  [N2](0,0)R(r)Z(z) 
is  assumed,  where  (r,z)  are,  respectively,  the  radial  and  axial  coordinates  and  there  is  no 
azimuthal  dependence.  Equation  (14)  can  then  be  written  in  the  form  showing  two  equations 
for  R  and  Z,  namely 


«odZ  D  d^Z  _  _D^  _  D_dR 
Z  dz  ^  Z  dz^  R  dr^  rR  dr 


(15) 


Since  the  two  equations  are  independent,  they  must  be  equal  to  a  common  constant  k„  where 
the  notation  for  the  species  in  level  v  is  reinstated.  Thus,  we  have 

%  sz  „ 

f  (16) 

dr^  r  dr  D 

The  appropriate  solutions  to  the  Eqs.  (16)  are  Z(z)  =  exp{k^)  and  R(r)  -Jq{F^  r/r^,  where 
Jq  is  the  zero  order  Bessel  function  of  the  first  kind  and  =  r^  Note  that  k^  and 

k^  are  wall  rates  in  this  development.  Substituting  the  solution  into  the  Z  equation  and  solving 
for  k„  gives 


^  _  _Uo_ 
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In  the  limit  of  k^  «  UqIAD,  Eq.  (17)  reduces  to  k„  =  -kju^.  Typically,  this  condition  was 
satisfied  to  better  than  two  orders  of  magnitude.  The  axial  fimction  now  becomes 
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Z(z)  =  exp(-A:^z/Mo)  =  exp(-^^0 


(18) 


where  t  =  zIuq  is  the  residence  time  of  the  gas  in  contact  with  the  specimen  wall.  Equation 
(18)  describes  the  exponential  decay  of  the  vibrationally-excited  nitrogen  gas  along  the 
specimen  tube  axis  due  to  deactivation  at  the  wall.  Substituting  the  solution  for  the  R 
equation  into  Eq.  (13),  gives  at  the  wall  with  r  =  Tq  the  expression 


(19) 


where  is  the  first  order  Bessel  function  of  the  first  kind.  If  the  deactivation  probability  y,, 
is  known,  then  Eq.  (1 9)  is  transcendental  in  wall  rate  (i.e.,  k^y,^.  However,  since  the  rate 
was  measured  directly  in  this  work,  Eq.  (19)  can  be  solved  for  y^  giving 


ADFyJ,{Fy) 

FrMF:)^2DFyJ,{Fy) 


(20) 


In  the  limit  of  small  y„ ,  Eq.  (20)  reduces  to  y„  =  For  example,  for  typical 

conditions  used  in  this  work,  for  a  value  of  =  1 5.5  s'\  this  approximation  is  8%  smaller 
than  the  value  of  y^  =5.07x10'^  given  by  Eq.  (20). 

Master  Equation  Model 

In  order  to  determine  the  impact  of  the  homogeneous  contributions  to  the  observed 
decay  rates,  a  well-tempered  theoretical  treatment  of  vibrational  energy  exchange  was 
explored.  Homogeneous  interactions,  specifically  the  V-V  processes,  mask  the  interactions 
with  the  wall.  The  heterogeneous  wall  loss  rates  determined  in  these  experiments  were 
greater  than  V-T  rates  by  about  three  orders  of  magnitude  while  the  V-V  rates  were  greater 
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than  the  wall  rates  by  a  similar  amount.  It  is  commonly  assumed  that  increases  linearly 
with  V.  Yet,  there  was  some  evidence  early  on  in  this  work  for  deviation  from  a  linear 
relationship  to  a  function  that  peaks  for  v  ^  4.  Moreover,  for  our  conditions,  the  V-V 
“pumping”  process  is  very  efficient  in  creating  a  nonBoltzmann  vibrational  population 
distribution  in  the  source  gas  that  is  subsequently  altered  by  the  interaction  with  the  wall. 
These  issues  indicate  that  there  will  be  a  difference  between  the  observed  decay  of  the 
vibrational  populations  (i.e.,  the  apparent  wall  rate)  and  the  actual  wall  rate.  Clearly,  as 
shown  in  Eq.(20),  the  deactivation  probability,  is  not  linearly  related  to  A:^„„,,)due  to  the 
diffusion  process.  Since  the  CARS  probe  is  located  in  the  gas  many  mean  free  paths  from  the 
surface,  it  can  be  expected  also  that  will  not  be  linearly  related  to  the  observed  decay  rate. 
For  these  reasons,  a  concerted  effort  was  made  to  fit  simulations  of  the  decay  of  vibrational 
populations  in  a  cylindrical  tube  based  on  the  Master  Equation  Model  to  all  deactivation  data 
in  order  to  account  for  the  impact  of  homogeneous  interactions  and  diffusion. 

The  simulation  code  was  previously  written  to  predict  the  vibrational  population 
distribution  of  nitrogen  under  only  the  influence  of  loss  terms  from  homogeneous  and 
heterogeneous  vibrational  modes.^  The  code  development  followed  previous  work  ofBailey^* 
and  Capitelli.^  Vibrational  energy  exchange  influences  certain  chemical  reactions  in  plasmas 
and  flames.  In  diatomic  molecules,  this  reactivity  increases  rapidly  with  v.  Theory^’  and 
experiment^*  show  the  dissociative  attachment  threshold  in  hydrogen  increases  by  five  orders 
of  magnitude  from  v  =  0  to  v  =  6.  The  dissociation  of  1%  silane  in  active  nitrogen  was 
measured  to  be  a  factor  of  70  larger  than  expected  from  theory.^’  Thus,  a  complete  picture 
of  the  dynamics  of  our  system  requires  more  detailed  knowledge  of  vibrational  pumping  and 
relaxation  processes.  The  V-V  exchange  rates  at  high  v's  are  particularly  difficult  to  measure 
because  the  rates  and  molecular  densities  are  extremely  small.  The  V-V  rates  in  the  lower  v's 
are  better  known  and  are  extrapolated  to  high  v's  in  the  simulation  code  by  the  well  known 
Schwartz,  Slawsky  and  Herzfeld  (SSH)  theory.'”’ 

The  V-V  exchange  rates  are  known  to  be  much  greater  than  the  V-T  rates  at  small  v: 
The  V-T  rate  increases  rapidly  with  v  and  surpasses  V-V  around  v  =  28  depending  on  the 
translational  temperature.  It  is  unlikely  that  the  V-T  interactions  had  any  direct  influence  on 


22 


the  population  distributions  observed  in  this  work,  but  they  were  included  in  the  code  to 
handle  the  upper  population  inversion  extending  to  v  =  47.  The  V-T  process  very  rapidly 
removes  population  for  v  >  47.  Translational  temperatures  in  the  measurements  were  always 
less  than  400  K.  These  temperatures  were  determined  by  recording  a  rotationally-resolved 
CARS  spectrum  for  the  ground  state  band  (v  =  0  to  v  =  1 )  during  each  experiment  and  fitting 
a  theoretical  CARS  spectrum  to  the  measured  spectrum.  The  theoretical  spectrum  was 
generated  by  the  Sandia  CARS  code.""  Rotational  populations  were  observed  to  be 
Boltzmann  distributed  in  all  vibrational  bands,  and  thus  were  equilibrated  to  the  translational 
temperature.  It  was  found  that  the  rotational  temperature  was  independent  of  v  consistent 
with  the  results  of  Massabieaux  et  al.‘*^  It  was  assumed  that  this  held  in  all  experiments 
carried  out  in  this  work.  With  the  modest  temperature  and  low  pressure  used  in  these  studies, 
three  body  collisions  rarely  occur.  Because  the  gas  temperature  was  somewhat  less  than  1 000 
K,  the  probabilities  for  multi-quantum  V-V  and  V-T  exchanges  were  very  small  compared  to 
single  quantum  steps,  and  thus  could  be  neglected."*^  All  heterogeneous  interactions  also  were 
calculated  on  a  single  quantum  loss  rate  basis.  The  model  assumed  the  molecule  to  be  a 
Morse-type  rigid  oscillator.  These  concepts  and  limitations  were  applied  to  the  Master 
Equation,  which  can  be  written  in  the  generalized  form 

Pi  =  S  (F^ ,  p*  -  F.  ^p.)  (21) 


here  r„„  is  the  conditional  probability  rate  of  the  transition  mton  and  P„  is  the  occupational 
probability  of  the  property  in  question.'^  Here  the  property  is  the  nitrogen  number  density, 
where  each  vibrational  level  is  taken  as  a  separate  species.  The  expression  used  in  the  model 
is  written^ 
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where  is  the  number  density  (cm'^)  of  a  neutral  nitrogen  molecule  in  the  \dbrational  state 
of  the  ground  electronic  state,  the  R^fs  are  the  exchange  rates  (s‘‘)  containing  the  V-T  and 
the  V-V  interactions  and  is  the  wall  deactivation  rate  of  the  v**  vibrational  state.  The 
rates  in  the  first  two  terms  on  the  right  side  of  Eq.  (22)  can  be  expanded  in  forms  given  by 
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which  represent  the  V-V  and  V-T  rates  of  single  quanta  energy  exchange  between  molecules 
wherein  population  is  introduced  into  level  v.  For  the  V-V  process,  population  increase  of 
the  V  level  given  in  the  second  term  of  the  first  expression  in  Eq.  (23)  derives  fi'om  one  level 
above  by  transferring  energy  to  a  second  molecule  and  the  corresponding  term  in  the  second 
expression  describes  the  population  increase  from  one  level  below  by  transfer  of  energy  from 
a  second  molecule,  where  the  are  the  rate  coefficients  (cmVs).  The  first  terms  on  the 
right  sides  of  the  two  expressions  in  Eq.(23)  correspond  to  V-T  processes  which  increase  the 
population  of  the  v  level,  where  the  P„„  are  the  rate  coefficients  (cmVs)  and  Vis  the  total 
number  density  of  the  gas.  The  corresponding  loss  rates  given  by  the  third  and  fourth  terms 
of  Eq.(22)  have  terms  analogous  to  Eq.(23). 
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ANALYSIS 


Data  Collection.  Corrections  and  Normalization 

Plots  of  typical  spectra  for  the  v  =  0,  1,  2  and  3  bands  are  shown  in  Fig.  3  The 
measurements  were  made  at  the  peaks  of  the  J  =  8  lines  of  the  vibrational  bands.  This 
simplified  both  the  deactivation  analyses  and  the  saturation  corrections;  however,  because  of 
the  narrow  width  of  the  rotational  lines,  the  CARS  signals  were  extremely  sensitive  to  the 
wavelength  setting  of  the  dye  laser.  Since  the  YAG  laser  was  operating  single  mode  (0.006 
cm'*),  the  observed  width  was  the  convolution  of  the  Doppler  lineshape  of  the  rotational  line 
with  the  lineshape  of  the  dye  laser.  Since  the  Doppler  width  was  0.006  cm'*  and  the  dye  laser 
lineAvidth  was  about  0. 1  cm'*,  the  observed  width  was  determined  by  the  dye  laser.  The  dye 
laser  wavelength  step  size  was  0.0005  nm  (0.0135  cm'*);  however,  the  resetability  was  two 
to  four  times  larger  than  this  step  size.  Thus,  it  was  very  difficult  to  dependably  move  from 
one  vibrational  band  to  the  next  without  misalignment  of  the  dye  laser  wavelength. 
Therefore,  the  rotational  line  in  each  vibrational  bands  was  scanned  and  then  fitted  to  the 
squared  Lorentzian  function.  This  procedure  provided  a  precise  measurement  of  the  peak 
CARS  intensity  with  a  much  relaxed  requirement  for  the  initial  wavelength  setting  of  the  dye 
laser.  This  type  spectral  measurement  is  termed  “skip  spectra.” 

Skip  spectra  were  obtained  in  both  deactivation  and  saturation  experiments.  In  a 
typical  experiment,  which  involves  several  runs  over  the  desired  vibrational  bands,  each  J  = 
8  line  is  scanned  with  16  spectral  points.  At  each  spectral  point  in  a  scan,  the  average  of  200 
to  500  laser  pulses  was  recorded.  For  each  laser  pulse,  the  PMT  produces  a  signal  CARS 
pulse  that  is  integrated  with  10-bit  precision  by  the  LeCroy  ADC.  The  signals  from  the  two 
laser-monitoring  PIN  diodes  are  also  integrated  with  the  ADC  and,  with  the  integrated  CARS 
signals,  are  averaged  at  each  spectral  point.  Signal  data  are  recorded  in  double  precision  for 
each  spectral  point  and  stored  in  files  as  sums.  This  includes  the  sum  and  the  sum  of  squares 
for  the  CARS,  dye  and  YAG  PIN  diode  signal  as  well  as  the  wavelength  and  number  of  pulses 
summed.  For  the  laser  covariance,  the  sums  of  sums  and  sums  of  differences  of  diode  signals 
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CARS  Signal  (a.u.)  CARS  Signal  (a.u.) 


and  their  sum  of  squares  are  also  written  to  the  file.  Separately,  hand  records  were  kept  of 
the  optical  attenuation  of  the  CARS  beam,  electronic  channel  gains  and  the  beam  powers 
between  each  run.  The  files  and  notes  were  used  to  determine  the  rescaled  spectral  value, 
saturation  correction  and  the  statistics. 

Rescaling  of  data  differs  between  a  saturation  and  a  deactivation  experiment.  For  the 
saturation  case,  each  band  v  is  normalized  separately  to  known  YAG  and  dye  laser  powers. 
In  the  deactivation  case,  an  entire  run  is  rescaled  using  the  average  PIN  diode  signals  from 
the  V  =  0  measurements  as  the  reference  values  and  then  the  raw  CARS  signals  are  corrected 
for  saturation.  The  rescaling  places  the  measurements  on  each  vibrational  band  in  a  given 
deactivation  experiment  on  the  same  the  basis  of  YAG  and  dye  powers  so  that  the  relative 
band  populations  can  be  extracted.  Examples  of  these  calculations  are 
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where  the  sums  are  for  each  spectral  point  i  of  the  band  v,  is  the  CARS  signal  sum,  82^^ 
is  the  sum  of  the  squares  of  the  CARS  signals,  Yj  is  the  sum  of  the  yag  PIN  diode  signals,  and 
Dj  is  the  sum  of  the  dye  PIN  diode  signals.  The  lower  case  quantities  are  the  pulse  values  and 
the  means  over  npts  pulses  are  calculated  later. 

Each  rotational  line  was  scanned  with  np  points  that  includes  two  points  in  the 
background.  The  CARS  background  signals  were  measured  close  to  the  band  heads  of  each 
vibrational  band  in  regions  essentially  fi'ee  of  any  detectable  rotational  lines.  The  background 
measurement  is  subtracted  from  the  average  signals,  which  removes  any  electronic  offset  and 
nonresonant  CARS  signal.  The  CARS  channel  offset  came  from  two  external  sources.  One 
was  a  slight  sensitivity  to  extraneous  electrical  emissions,  principally  from  the  flash  lamp 
power  supply,  and  the  other  was  leakage  of  the  532  ran  light  onto  the  PMT.  Electronic 
offsets  in  PIN  diode  channels  did  not  change  more  than  five  counts  in  three  years  and  was. 
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thus  constant  throughout  an  experiment.  Offsets  in  the  PIN  diode  channels,  and  , 
were  determined  by  observing  the  dark  count  value  for  each  detector,  which  was  typically  25 
counts. 

CARS  background  exhibited  a  slight  gradual  change  across  a  vibrational  band.  This 
characteristic  was  constant  through  an  experiment  and  probably  over  many  experiments. 
Actually,  the  nonresonant  background  was  weak  or  nonexistent.  A  conventional  background 
correction  uses  measurements  on  either  side  of  the  feature  to  predict  the  background  at  the 
peak  of  the  line.  In  this  work,  the  value  on  one  side  was  subtracted,  which  provided  a 
consistent  and  expedient  method.  The  error  was  usually  well  under  10  counts  or  1%.  This 
error  was  largest  in  the  high  hot  bands  where  the  signal-to-background  ratio  was  the  lowest. 


The  band  averages  for  the  dye  laser  and  YAG  laser  PIN  diode  signals  over  np  spectral 
points  are  computed  from 
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where  and  Y^  were  25,  and  the  rescaled  CARS  signal,  ,  was  calculated  for  each 
spectral  point  i  in  band  v  as  given  in  by 
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which  normalizes  each  spectral  datum  to  the  YAG  and  dye  power  averages  for  that  scan.  The 
variances  of  the  means  of  the  CARS  signals,  the  dye  and  YAG  PIN  diode  signals  and  the 
covariance  between  the  dye  and  the  YAG  diode  signals  were  also  calculated.  The  correlation 
coefficient  between  diode  signals  was  typically  between  0.2  and  0.1  or  less  which  supports 
the  normalization  using  the  averages. 

Saturation  corrections,  derived  from  calibration  experiments,*^  were  used  to  correct 
the  reduction  in  response  caused  by  stimulated  Raman  scattering  (SRS)  which  pumps 
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population  directly  into  the  upper  level  simultaneously  vvith  the  CARS  process.'*’ 
Saturation  calibration  measurements,  which  gave  the  relationship  between  the  peak  CARS 
signals  and  laser  powers,  would  be  conducted  typically  within  no  more  than  20  days  of  the 
deactivation  experiments  to  which  the  calibration  was  applied  providing  that  the  optical 
alignment  of  the  system  had  not  been  disturbed.  The  average  dye  laser  power  was  varied  in 
steps  from  ~33  to  ~  10  mW  for  the  saturation  measurements.  Usually,  30  mW  (1.5  mJ/pulse 
at  20  Hz)  was  the  power  used  in  the  deactivation  experiments.  Power  incident  on  the  main 
chamber  was  a  convenient  unit  for  the  calibration  of  two  PIN  diodes.  The  J  =  8  line  was 
scanned  in  at  least  three  of  the  normally  detected  hot  bands  of  v  =  1, 2, 3,  and  4  at  each  dye 
laser  power  step  while  holding  all  other  system  variables  constant.  Nine  to  12  power  steps 
were  recorded.  Using  the  procedure  described  by  Yaney  and  Parish,'*  the  average  CARS 
signal  (PMT  counts)  in  band  v,  at  the  J  =  8  peak  and  the  laser  powers  and  (PIN 
diode  counts)  were  combined  to  plot  (So2v/Dv  ’Yy^)"*  versus  ’Y^  for  each  v-band. 
P^ameterized  in  this  way,  the  data  gives  a  linear  relationship  for  each  vibrational  band  with 
an  intercept  b,,  at  zero  laser  power  and  a  slope  m^  as  given  in 
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where  m^  is  negative.  The  square  root  of  the  ratio  of  the  unsaturated  to  saturated  signal,  8,3,  v, 
gives  the  correction  factor  for  each  vibrational  level  as'* 
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The  ratio  m^bv  in  Eq.(28)  generally  is  a  linear  function  of  v.'*  Therefore,  doing  a  linear 
regression  of  these  ratios  to  v  +  1  from  as  few  as  three  v  bands  provided  a  convenient  way 
to  accomplish  the  correction  of  the  deactivation  data  for  all  v  bands  under  each  particular  set 
of  experimental  circumstances.  Since  the  PIN  diode  responses  to  the  powers  of  their 
respective  lasers  had  been  calibrated,  the  saturation  calibrations  could  be  used  to  correct  the 
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deactivation  data  using  count  values  in  Eq.(28).  It  can  be  shown  that  the  methods  used  are 
equivalent  to  using  power  values.  The  deactivation  and  saturation  experiments  were 
connected  by  diode  counts-per-mW  (cnts/mW)  ratios  which  were  used  to  convert 
deactivation  power  readings  to  appropriate  saturation  values. 

The  saturation  measurements  were  normally  carried  out  by  varying  the  dye  laser 
power  by  rotating  a  halfwave  plate  just  in  front  of  a  polarizer  while  holding  the  YAG  power 
constant.  To  check  that  the  saturation  effect  depended  on  the  product  of  the  two  laser 
intensities  in  the  probe  volume  as  shown  in  Eq.(28),  the  YAG  power  was  varied  in  one 
experiment  while  holding  the  dye  power  constant.  The  CARS  signal  had  the  same  saturation 
dependence  on  the  incident  YAG  power  as  it  did  by  varying  the  dye  power. 

The  data  at  the  spectral  points,  are  free  of  offset,  background  and  any  laser  power 
drift  during  the  scan  over  the  line.  These  spectral  data  of  the  rotational  lines  were  fitted  to 
a  Lorentzian  squared  fianction  to  obtain  the  peak  heights  as  discussed  in  the  next  section.  The 
rotational  peak  heights  were  then  corrected  for  saturation  and  the  resulting  values  were 
rescaled  to  the  average  YAG  and  dye  laser  powers  present  during  the  measurements  on  the 
V  =  0  to  1  ground-state  band  in  that  experiment.  However,  the  optical  attenuation  and  the 
amplifier  gain  were  adjusted  for  each  hot  band  to  fit  the  observed  signals  well  within  the 
1 024-count  digitization  range  of  the  ADC.  The  dynamic  range  of  signals  over  the  four  or  five 
vibrational  bands  was  greater  than  7x10*.  This  dynamic  range  was  covered  by  the  neutral 
density  filter  wheel  that  provides  attenuations  from  unity  to  2145  and  the  amplifier  that  has 
10  gain  steps  from  unity  to  10.  Generally,  the  attenuation  and  the  gain  were  adjusted  to  keep 
the  maximum  CARS  signal  pulses  within  750  counts  or  less  of  the  ADC.  The  CARS  signals 
displayed  considerable  fluctuations  attributed  to  the  mode  beating  in  the  dye  laser  pulse  which 
made  these  adjustments  very  critical  to  acquiring  viable  data. 

Fitting  the  rotational  line 

In  order  to  achieve  a  high  degree  of  precision  in  determining  the  CARS  sign^ 
strength,  the  J  =  8  rotational  line  in  each  vibrational  band  was  measured  by  scanning  this  line 
using,  typically,  14  steps  of 0.0005  ran  (0.0135  cm'*)  and  200  laser  pulses  per  step  plus  two 
background  steps.  The  background  signal  level  was  removed  from  these  raw  data  and  the 
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square  of  the  Lorentzian  profile  was  used  to  fit  to  the  background-free  data  using  nonlinear 
regression  analysis.  This  profile  can  be  written  for  narrow  lines  in  wavelength  form  as 


(29) 


where  Wj  is  the  square  root  of  the  peak  CARS  amplitude  of  the  line,  Wj  is  the  wavelength  at 
the  peak,  Wj  is  the  linewidth  (FWHM)  and  A,;  is  the  wavelength  at  the  i*  point  of  the  scan. 
Because  any  background  signal  was  removed  prior  to  applying  Eq.  (29),  no  additive  constant 
was  needed  in  the  expression.  Parameter  Wj  is  proportional  to  the  population  difference 
N(v,  J  =  8)  -  N(v  + 1,  J  =  8),  which  was  then  used  to  compute  the  relative  populations  among 
the  observed  vibrational  bands.  The  fitted  values  for  Wj  were  recorded  in  a  disk  file  together 
with  the  PIN  diode  count  averages,  and  .  These  values  of  w,  were  then  fiilly  corrected, 
as  stated  at  the  end  of  the  last  section,  to  the  ground  state  laser  powers  and  for  saturation. 
The  corrected  values,  w^,,  for  band  v  are  obtained  from 
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where  Yq  and  Dq  correspond  to  the  average  ground  state  (v  =  0  to  v  =  1)  laser  power  values, 
G  is  the  amplifier  gain  and  A  is  the  optical  attenuation  used  in  the  band  v  measurements. 

Obtaining  the  observed  deactivation  rates 

Following  the  various  normalizations  and  corrections  of  the  peak  values  of  the  J  =  8 
line  in  each  vibrational  band,  the  data  for  the  v  bands  at  each  residence  time  were  processed 
together  to  extract  fractional  populations.  The  residence  time  is  simply  the  length  of  the 
specimen  tube  exposed  to  the  vibrationally-excited  gas  divided  by  the  gas  flow  velocity  in  the 
tube.  It  is  important  to  emphasize  that  the  peak  line  strengths,  w^,, ,  obtained  at  a  given 
residence  time  were  processed  together  independent  of  the  data  at  other  residence  times.  As 
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shown  in  Eq.(8),  the  line  strength  is  proportional  to  the  CARS  susceptibility,  Xcars,  which 
is  given  by'* 
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where  the  frequency  Ogj,  is  the  frequency  of  the  transition  between  levels  a  and  b,  daldQ.  is  the 
differential  cross  section  of  the  linear  Raman  process  which  depends  on  v,  p  and  s  identify 
pump  and  probe  (Stokes)  quantities,  h  is  Planck’s  constant  divided  by  27C,  and  Nb  are  the 
population  densities  and  F^b  is  the  line  width.  The  total  susceptibility  includes  the  real  and 
imaginary  parts  obtained  from  Eq.(31)  and  a  nonresonant  component  as  shown  in  Eq.(32). 

Xfor  X„r  (32) 

If  the  nonresonant  component  has  sufficient  magnitude  it  can  complicate  the  functional  form 
of  the  CARS  spectrum  due  to  the  squaring  in  Eq.(8).  In  this  work,  Xm  was  essentially  zero. 
At  resonance  when  Wp  -  o),  =  tDjb»  which  corresponds  to  the  peak  of  the  line,  only  the 
imaginary  part  of  Eq.(3 1)  remains.  Thus,  using  the  peak  line  strength  simplifies  the  analysis. 
It  is  assumed  in  this  analysis  that  the  rotational  temperature,  and  therefore  the  rotational 
population  distribution,  is  equilibrated  with  the  translational  or  gas  temperature.  The 
population  distribution  over  v  and  J  is  adequately  described  by 

(33) 

where  N  is  the  total  density  of  molecular  nitrogen,  Z  is  the  partition  function,  G(v)  is  the 
vibrational  term  energy  from  the  ground  state,  v  and  J  are,  respectively,  the  vibrational  and 
rotational  quantum  numbers,  Bg  is  the  equilibrium  rotational  constant,  is  the  vibration- 
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rotation  interaction  constant,  T„  is  the  vibrational  temperature  of  band  v  corresponding  to  the 
V  to  V  + 1  transition,  is  the  rotational  temperature,  k  is  the  Boltzmann  constant  and  gj  is  the 
statistical  weight  due  to  nuclear  spin  degeneracy  which  depends  on  J  for  nonzero  nuclear  spin. 
The  partition  function  is  found  by  summing  Eq.(33)  over  v  and  J;  however,  it  can  be  factored 
into  vibrational  and  rotational  parts  with,  at  most,  a  1%  error  by  neglecting  the  right  hand 
factor  involving  v  and  J  in  Eq.(33).  Thus, 
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and  Zf  is  the  sum  of  the  J-only-dependent  factors  in  Eq.(33).  The  population  difference  for 
a  Q  Branch  transition  where  AJ  =  0,  AN.^,  j  =  N„  j  -  Nv+,,j  can  be  written  as 

where  K^j  =  exp[a^(v+0.5)J(J+l)/A7’J 
and  K'j  =  ^;(2y+ l)exp[-fi/(J+ 1)/A7;]/Z^ 


From  Eqs.(8)  and  (33),  the  square  root  of  the  CARS  signal  can  be  written 
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where  is  the  Raman  cross  section  of  the  v  vibrational  band  and  h(v)  corrects  for  the 
deviation  of  the  Morse  potential  from  an  harmonic  oscillator.^  The  width  of  the  Doppler- 
broadened  rotational-vibrational  line  depends  on  the  frequency  of  the  transition  and  the  gas 
temperature.  This  introduces  another  correction  factor,  ,  given  by'** 
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The  strength  of  the  transition  is  inversely  proportional  to  the  width.  Ey  j  is  the  energy  of  the 
rotational- vibrational  level  and  Eq  .(3  7)  references  the  strength  of  a  transition  in  band  v  to  the 
ground  state  band  v  =  0.  All  measurements  were  made  on  J  =  8.  Combining  Eqs.(35)-(37) 
assuming  constant  laser  powers,  perfect  phase  matching  and  constant  beam  intersection 
geometry,  the  square  root  of  the  CARS  signal  becomes 

'  C,(v  t  JV„,)  (38) 


where  Cj  contains  the  laser  power  factors  and  all  the  constants  in  Eqs.  (8),  (31),  and  (34)- 
(36).  At  some  v  band  the  population  will  be  so  small  that  the  signal  will  be  unmeasurable. 
It  is  assumed  that  populations  at  higher  v’s  had  negligible  affect  on  the  measured  populations 
in  the  lower  v  bands.  This  permitted  a  truncation  of  the  calculation  for  extraction  of  fractional 
populations.  Thus,  by  solving  Eq.(38)  for  Ny  gives  a  summation  from  v  to  the  maximum  v 
measured,  v^ ,  which  for  J  =  8,  )delds  Eq.(39),  where  Ny’  =  NyCj  is  a  relative  population. 
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In  this  summation,  the  population  of  the  Vn,,^  +1  band  is  set  to  zero.  The  total  relative 
population,  N’,  is  then  the  sum  of  Eq.(39)  from  v  =  0  to  v  =  v„„.  Thus,  the  fractional 
populations,  Fy,  are  found  from 
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Various  schemes  were  explored  to  extrapolate  the  calculation  of  Fy  to  include  higher  values 
of  v^  based  on  theoretical  projections.  Although  these  extrapolations  produced  small 
changes  in  the  values  of  F4  and  F5  compared  to  the  above  truncation  procedure,  these  changes 
had,  at  most,  only  a  small  impact  (<10%)  on  the  observed  deactivation  rates. 
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DATA  AND  RESULTS 


Tables  1-6  summarize  the  results  of  this  project.  Reproducibility  of  these  results 
varied  because  of  accumulated  surface  changes  due  to  use  of  the  specimen  in  many 
experiments,  planned  treatments  of  the  specimens,  undetermined  variations  in  the  specimens 
or  uncontrollable  variations  in  the  experimental  system.  Several  experiments  with  P)n'ex 
tested  these  possibilities.  Tables  1,  3  and  5  report  the  observed  deactivation  rates,  ,  and 
Tables  2, 4  and  6  give  the  corresponding  apparent  deactivation  coefficients,  computed 
directly  from  the  (viz.,  instead  of  k^^,)  values  using  Eq.(20).  Generally,  is  expected 
to  increase  for  higher  pressure  and  decrease  with  increased  temperature.  For  a  given 
pressure,  the  temperature  for  specified  discharge  conditions  was  determmed  by  the  gas  flow 
rate.  Higher  flow  rates  gave  higher  temperatures  due  to  increased  heat  transport.  The 
observed  deactivation  rate,  however,  is  a  combination  of  the  interactions  of  the  molecules 
with  the  wall  and  the  interactions  with  other  molecules  through  the  V-V  exchange  process. 
This  greatly  complicates  the  interpretation  of  the  data.  The  fitting  of  the  Master  Equation 
Model  attempts  to  extract  the  wall  deactivation  rate,  but  the  interplay  between  these  two 
interactions  can  vary  significantly  from  experiment  to  experiment.  The  results  of  these  fittings 
are  discussed  in  the  next  section. 

Except  for  gold,  all  surfaces  were  found  to  be  nitrided  by  the  afterglow  due  to  the 
small  density  of  nitrogen  atoms.*^  This  density  was  estimated  to  be  0.55%  of  [N2]  from 
measurements  of  the  decay  of  the  first  positive  emission,  B^IIg  to  A^S„,  versus  residence  time. 
Surface  measurements  on  Pyrex  using  XPS  showed  nitrogen  bonded  principally  to  the  silicon 
and  secondarily  to  the  oxygen  surface  atoms  of  witness  coupons.^  The  complete  XPS  results 
are  reported  in  a  later  section.  Coverage  was  determined  to  be  between  8%  and  10%  of  the 
surface  with  one  case  showing  about  30%.  It  is  clear  that  exposure  to  the  afterglow  changes 
the  surface  in  some  way.  The  details  of  the  bonding  mechanisms  cannot  be  obtained  nor 
completely  understood  from  these  experiments.  Nitrogen  adatoms  are  expected  to  make  the 
surface  passive  by  permanently  occupying  surface  adsorption  wells  generated  by  dangling 
bonds.  They  also  change  the  principal  surface  resonance  vibration  away  from  that  of  the 
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virgin  surface  which  could  provide  more  efficient  coupling  to  vibrating  molecules.  Thus, 
adatoms  may  either  increase  or  reduce  the  surface  deactivation  rate  and,  in  the  same  respect, 
the  number  of  surface  collisions  a  molecule  makes  or  the  number  of  quanta  lost  per  collision 
may  be  impacted  by  the  nitriding. 


Pvrex 

Table  1  gives  the  observed  values  of  k^,  with  fitting  uncertainties  and  Table  2  gives  the 
corresponding  apparent  deactivation  coefficients  Yva-  These  tables  also  contain  results  of 
Master  Equation  modeling  efforts.  In  Table  1,  “nitrided”  indicates  the  surface  had  been 
exposed  to  active  nitrogen  (i.e.,  afterglow)  in  several  previous  experiments  and  was  found 
from  XPS  measurements  to  have  nitrogen  bonded  to  the  surface;  “neat”  refers  to  the  first  time 
this  specimen  was  used,  taken  as  is  from  the  glass  shop  and  treated  only  to  a  general  cleaning; 
“+column”  refers  to  an  experiment  where  the  tube  was  placed  in  the  positive  column  for  33.7 
hours  before  the  deactivation  measurements;  and  “exposed”  refers  to  an  experiment  where 
the  tube  spent  39.5  hours  in  the  afterglow  before  deactivation  measurements.  The  cleaning 
procedure  generally  consisted  of  ultrasonic  degreasing  in  methylene  chloride  followed  by  ion 
cleaning  at  about  5  Torr  in  flowing  argon  for  30  to  60  minutes. 

Experiments  1  and  2  used  the  same  tube  wherein  the  latter  experiment  occurred  after 
additional  exposure  of  the  tube  to  the  afterglow  thereby  increasing  the  degree  of  nitriding. 
It  would  seem  that  the  additional  nitriding  in  Experiment  2  is  responsible  for  the  higher  values 
of  Yva  The  lower  flow  rate  and  the  resulting  lower  temperature  introduces  ambiguity  in  any 
attempt  to  interpret  these  data.  Changes  due  to  exposure  to  the  afterglow  were  also  reported 
by  Black  et  al.  using  a  microwave  source*;  however,  they  report  a  decrease  of  Yi  for  Pyrex 
after  a  24  hour  exposure  to  the  afterglow.  They  found  no  change  with  exposure  when  using 
a  thermal  source.  This  can  be  explained  by  the  fact  that  they  measured  only  the  v  =  1  band 
using  a  1-m  long  specimen  tube.  With  the  microwave  source,  as  with  the  positive  column 
source  used  in  this  work,  the  vibrational  population  is  distributed  across  all  low-v  states 
corresponding  to  a  vibrational  temperature  typically  well  above  the  2000  K  they  report  for 
their  thermal  source.  Moreover,  there  is  no  surface  nitriding  using  a  thermal  source.  Thus, 
measuring  only  the  v  =  1  population  produced  by  a  discharge  source  is  not  sufficient  to 
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correctly  characterize  the  population  changes  due  to  wall  deactivation  over  long  residence 
times,  especially  with  V-V  pumping  and  nitriding  effects. 

Experiments  3  and  4  were  carried  out  to  check  for  reproducibility  and  the  degree  of 
dependence  on  pressure.  The  data  and  fitting  results  for  Experiment  3  are  given  in  Fig.  4. 
The  flow  rate  in  Experiment  4  was  not  increased  to  compensate  for  the  higher  pressure; 
therefore,  the  temperature  in  Experiment  4  was  not  the  same  as  in  Experiment  3. 
Nevertheless,  the  deactivation  rates  in  these  two  experiments  are  the  same  within  fitting 
uncertainties.  Experiment  5  attempted  to  achieve  a  high  degree  of  nitriding.  Although  the 
Yyj  values  are  larger  than  the  values  obtained  for  the  “neat”  specimens,  the  rate  of  increase 
of  Yva  with  V  for  the  “+  column”  specimen  is  about  half  that  of  the  neat  specimens. 
Experiments  4,  6  and  7  used  one  tube  and  Experiment  7  followed  the  next  day  after 
Experiment  6.  Experiments  6  and  7  were  run  to  test  reproducibility.  The  Yva  values  show 
reasonable  agreement  within  the  fitting  uncertainties  plus  the  estimated  5  to  10%  uncertainty 
in  the  repeatability  of  the  set  up  of  the  experimental  conditions.  Comparison  of  Experiments 
6  and  7  to  Experiment  4  again  shows  the  rate  of  increase  of  Yva  with  v  to  be  slower  after  the 
tube  was  exposed  to  the  afterglow.  An  experiment  was  attempted  at  a  pressure  of  0.5  Torr 
also  on  the  this  specimen.  Although  the  signals  were  very  weak,  the  observed  deactivation 
rates  for  v  =  1  and  2  were  about  half  the  values  obtained  in  Experiments  6  and  7  confirming 
the  expected  pressure  dependence. 


Gold  Foil 

The  results  reported  in  Tables  2  and  3  are  for  the  same  gold  foil  and  are  tabulated  in 
order  of  acquisition.  The  foil  was  formed  to  completely  cover  the  inside  of  a  Pyrex  tube. 
Gold  was  chosen  because  surface  changes  due  to  reaction  with  active  nitrogen  were  not 
expected.  Thus,  experiments  could  be  conducted  to  assess  the  effects  of  temperature 
(determined  by  the  gas  flow  rate)  and  pressure  in  the  absence  of  nitriding.  It  should  be  noted, 
however,  that  the  surface  did  change  from  a  calendered  gloss  to  a  slightly  darker,  very  fine 
buff.  This  surface  roughening  could  only  be  caused  by  collisions  of  vibrationally  hot  nitrogen 
during  the  experiments  or  by  the  argon  positive  column  used  for  surface  cleaning.  The  first 
three  experiments  in  Tables  3  and  4  used  three  different  gas  flow  rates  in  order  to  vary  the  gas 
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Fractional  Population 


Figure  4.  Deactivation  measurements  on  neat  Pyrex  at  360  K  and  1 1  Torr.  The  points  are  the  data; 
the  solid  lines  are  expontential  fits  to  the  data  with  the  rates  given;  the  dashed  lines 
are  the  results  of  the  fitting  of  the  Master  Equation  Model  to  the  data;  and  the 
vibrational  temperature  of  the  v  =  1  band,  Ti,  is  given  by  the  triangles  joined  by  the 
cubic  spline  curve. 
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temperature.  The  tube  was  mounted  in  the  chamber  using  glass  beads  that  provided  good 
thermal  isolation.  Thus,  it  is  assumed  that  the  tube  temperature  was  equal  to  the  gas 
temperature  as  measured  from  a  CARS  rotational  spectrum.  As  expected,  the  deactivation 
rates  decreased  uniformly  with  increase  in  temperature.  The  last  three  experiments  were 
intended  for  pressure  studies.  Comparing  Experiments  1  and  4,  which  are  at  the  same 
temperature,  show  significant  decreases  in  deactivation  rates  and  the  corresponding  Yva  values 
with  drop  in  pressure;  however,  the  opposite  change  is  shown  in  Experiments  2  and  5  which 
were  performed  at  a  higher  temperature.  It  is  likely  that  this  apparent  contradiction  is  a 
consequence  of  the  complexities  of  the  mixing  of  the  wall  deactivation  effects  and  the  V-V 
interactions.  The  data  and  fitting  results  for  Experiment  2  are  given  in  Fig.  5. 

Selected  Materials 

Tables  5  and  6  report  the  results  on  a  variety  of  materials  in  the  same  form  as  the 
previous  tables.  The  two  stainless  steel  specimens  showed  the  highest  deactivation  rates  as 
would  be  expected  and  as  reported  by  Black  et  al.*  The  304  alloy  was  measured  at  the  start 
of  these  studies  when  the  discharge  tube  was  operated  at  a  higher  temperature;  however,  the 
inside  surface  of  this  tube  had  been  machined  to  a  fairly  rough  surface.  On  the  other  hand, 
the  AMS  312  tube  had  a  very  smooth  “drawn”  surface.  Although  the  alloys  are  different,  the 
high  deactivation  rates  obtained  with  the  304  specimen  is  attributed  to  the  large  effective 
surface  area.  Figure  6  shows  plots  of  the  data  and  the  fitting  results  for  AMS  312.  Two  sets 
of  values  are  given  for  aluminum  which  represent  experiments  on  two  different  tubes  of  the 
same  alloy.  The  results  are  in  reasonable  agreement  which,  again,  shows  that  the  measurement 
procedures  were  carried  out  in  a  consistent  manner.  The  results  for  titanium  were  a  surprise. 
The  same  tube  was  used  in  both  experiments.  Experiment  3  was  first  and  showed  virtually  no 
deactivation  of  the  v  =  1  band  while  Experiment  4  showed  that  the  surface  had  been  altered 
by  the  exposure  to  active  nitrogen.  This  nitriding  again  reduced  the  rate  of  increase  of 
deactivation  rate  with  v,  which  is  a  factor  of  two  between  Experiments  3  and  4 .  Titanium  gave 
the  lowest  deactivation  rates  observed  in  these  studies,  which  has  to  be  attributed  to  the  oxide 
layer  on  the  titanium  surface.  Teflon,  chosen  to  represent  polymer-based  materials,  was 
expected  to  be  benign.  This  material,  indeed,  gave  the  lowest  deactivation  rates 
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Figure  5.  Deactivation  measurements  on  gold  foil  at  351  K  and  17  Torr.  The  points  are  the  data; 
the  solid  lines  are  expontential  fits  to  the  data  with  the  rates  given;  the  dashed  lines 
are  the  results  of  the  fitting  of  the  Master  Equation  Model  to  the  data;  and  the 
vibrational  temperature  of  the  v  =  1  band,  Ti,  is  given  by  the  triangles  joined  by  the 
cubic  spline  curve. 
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Figure  6.  Deactivation  measurements  on  AMS  312  stainless  steel  at  340  K  and  14  Torr.  The 
points  are  the  data;  the  solid  lines  are  exponential  fits  to  the  data  with  the  rates  given; 
the  dashed  lines  are  the  results  of  the  fitting  of  the  Master  Equation  Model  to  the 
data;  and  the  vibrational  temperature  of  the  v  =  1  band,  Ti,  is  given  by  the  triangles 
joined  by  the  cubic  spline  curve. 
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except  for  titanium.  With  two  fluorine  atoms  per  carbon  atom.  Teflon  contains  a  high  fluorine 
fraction  which  would  otherwise  contribute  to  strong  deactivation.  The  stability  of  these  bonds 
is  generally  quite  high,  and  therefore  Teflon  would  be  expected  to  have  only  a  weak  interaction 
with  vibrational-excited  species.  Nevertheless,  it  was  found  very  important  to  keep  the  gas 
temperature  low  with  this  material,  otherwise  the  surface  would  be  altered  and  release 
contaminates  into  the  flow. 


Table  1:  Observed  deactivation  rates  k^  for  Pyrex 


Exp 

Surface 

T 

(K) 

P 

(Torr) 

Gas  Flow 
(seem) 

^1 

(s-‘) 

^3 

(s-‘) 

Ml 

1 

nitrided 

315 

17 

600 

7.2 

±6% 

20.2 

±4% 

24 

±7% 

26 

±10% 

2 

nitrided 

307 

17 

450 

12.9 

±3% 

28 

±3% 

37 

±7% 

n 

3 

neat 

360 

11 

802 

12 

±14% 

30 

±7% 

47 

±5% 

58 

±9% 

■ 

neat 

330 

14 

800 

11 

±12% 

31 

±5% 

46 

±2% 

62 

±5% 

5 

+column 

321 

11 

800 

19 

±11% 

38 

±5% 

57 

±5% 

74 

±11% 

6 

exposed 

325 

14 

800 

14.9 

±5% 

34 

±5% 

47 

±4% 

61 

±4% 

■ 

exposed 

328 

14 

800 

14 

±7% 

^9 

42 

±7% 

51 

±9% 

Table  2.  Apparent  deactivation  coefficients  Yv,  1 0^  for  Pyrex 


Exp 

1  Yu 

y2a 

Y3a 

Y4a 

1 

9.5 

10 

2 

4.6 

11.8 

18 

1 

1 

3 

9.7 

17.4 

21 

n 

3.5 

11.9 

19.7 

32 

5 

W 

17.2 

28.3 

35 

6 

5.1 

13.6 

21 

32 

■ 

11 

18 

24 

The  %uncertainty  in  «  (l+0.01kv)x%uncertainty  in  k^ 
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Table  3:  Observed  deactivation  rates  of  gold  foil 


Exp 

T 

(K) 

P 

(Torr) 

Gas  Flow 
(seem) 

ki 

(s-') 

kj 

(s’*) 

ks 

(s-‘) 

k4 

(s-*) 

ks 

(s-') 

1 

323 

17 

600 

22 

±4% 

51 

±2% 

62 

±3% 

73 

±3% 

-  -  - 

2 

351 

17 

900 

17 

±5% 

37 

±3% 

53 

±3% 

64 

±3% 

65 

±5% 

3 

379 

17 

1200 

7.5 

±9% 

19 

±6% 

32 

±4% 

42 

±4% 

50 

±4% 

■ 

323 

10 

900 

10 

±11% 

20 

±4% 

34 

±5% 

42 

±8% 

48 

±6% 

5 

353 

23 

900 

15.7 

±6% 

28 

±3% 

48 

±6% 

60 

±9% 

— 

6 

340 

14 

800 

! 

20.7 

±4% 

47 

±7% 

67 

±4% 

83 

±4% 

Table  4.  Apparent  deactivation  coefficients  x  10'*  of  gold* 


Exp 

Yu 

y2a 

Ysa 

Y4a 

Ysa 

1 

8.2 

25.9 

38 

57 

— 

6.6 

14.9 

25 

33 

32 

2.2 

6.0 

11.0 

15.4 

20 

D 

3.1 

6.8 

12.1 

16 

19 

5.3 

11.0 

24 

41 

-  -  - 

7.0 

19 

33 

49 

— 

*  The  %uncertainty  in  y„  «  (l+0.01ky)x%uncertainty  in 
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Table  5.  Observed  deactivation  rates  of  selected  materials 


Exp 

Material 

T 

(K) 

P 

(Torr) 

Gas  Flow 
(seem) 

H 

k2 

(s-‘) 

ks 

(s-*) 

M 

1 

Stainless  steel 

460 

17 

1000 

53.1 

108 

138 

S304 

B 

2 

Stainless  steel 

340 

14 

800 

19 

44 

70 

135 

AMS  312 

±9% 

±5% 

±6% 

±7% 

3 

Aluminum 

392 

17 

1000 

12 

a 

31 

38 

3003 

349 

14 

800 

8.3 

HI 

32 

42 

±7% 

±7% 

±6% 

±6% 

D 

Titanium 

382 

14 

1100 

«0 

5 

8 

13 

II 

AMS  4943D 

±26% 

±21% 

±20% 

Titanium 

375 

14 

800 

1.3 

3.7 

6.4 

8.1 

AMS  4943D 

±22% 

±11% 

±9% 

±8% 

6 

Teflon 

334 

14 

800 

3.4 

8.8 

14 

20 

±15% 

±10% 

±7% 

±8% 

Table  6.  Apparent  deactivation  coefficients  Yva  ^  10**  of  selected  materials 


Exp 

Surface 

Yla 

Yaa 

Yaa 

Y4a 

1 

SS  304 

18 

53 

96 

— 

2 

SS312 

7.8 

23 

58 

10" 

3 

A13003 

4 

7.0 

wm 

15 

2.8 

7.3 

Eg 

18.3 

a 

Ti  4943D 

«0 

1.5 

2.5 

4.3 

5 

Ti  4943D 

0.4 

1.2 

2.1 

2.9 

6 

Teflon 

1.1 

2.8 

wm 

7 

The  %uncertainty  in  yv,  *  (l+0.011cv)x%uncertainty  in  k,, 


Master  Equation  Modeling 

The  values  given  in  Tables  2,  4  and  6  include  the  influence  of  the  V-V  processes  as  well  as 
interactions  with  the  wall.  The  Master  Equation  Model  (MEM)  computes  the  population 
distributions  up  to  v  =  40;  however,  the  measurements  extend  no  higher  than  v  =  5  and  more 
often  only  to  v  =  4.  In  general,  above  some  value  of  v,  the  deactivation  coefficient  is  unity, 
which  occurs  from  as  low  as  v  =  4  to,  perhaps,  v  =  10.  Thus,  it  was  necessary  to  devise  a 
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reasonable  procedure  for  extrapolating  the  to  higher  v’s  in  order  to  use  the  MEM  program. 
Two  methods  were  used:  1)  since  the  observed  deactivation  rates  typically  increased  linearly 

with  V,  a  linear  extrapolation  of  k^  was  carried  out  until  reached  unity;  2)  the  square  root 
of  the  measured  CARS  signals  were  extrapolated  using  the  Treanor  distribution  concatenated 
on  to  the  experimental  curve  using  the  two  highest  v-band  populations  measured  at  each 
residence  time.  The  latter  approach  creates  a  renormalization  of  the  data  such  that  the 
recorded  values  of  1^,  are  shifted  slightly  from  the  values  given  in  Tables  1, 3  or  5.  In  either 
approach,  the  actual  wall  deactivation  rates  k^vvaii)  to  the  data  from  which  the  “true” 

deactivation  coefficients  y^  are  calculated  using  Eq.(20).  Table  7  gives  the  results  of  six 
studies  of  fitting  the  MEM  to  the  data.  It  should  be  noted  that  these  calculations  were  very 
tedious  and  the  fittings  were  often  plagued  by  ambiguities.  Because  of  the  shear  complexity 
of  the  relationships  between  the  deactivation  rates  of  the  observed  vibrational  states  and  the 
impact  of  the  interplay  between  the  V-V  and  wall  deactivation  rates  on  the  radial  distribution, 
it  was  often  impossible  to  find  a  unique  and  plausible  fit  to  a  given  data  set. 

The  results  in  Table  7  show  a  very  significant  characteristic.  A  surface  giving  high 
wall  deactivation  rates  is  accompanied  by  observed  deactivation  rates  that  are  lower  than  the 
wall  values.  This  characteristic  also  carries  over  to  the  deactivation  coefficients  wherein  the 
Yv  are  larger  than  the  apparent  values  y^.  This  can  be  seen  for  the  +  column  Pyrex,  the 
stainless  steel  and  the  aluminum  results  given  in  Table  7.  On  the  other  hand,  the  materials 
giving  low  or  weak  wall  effects  show  wall  deactivation  rates  I^^id  that  are  very  similar  to  or 
smaller  than  the  observed  deactivation  rates  ky.  Likewise,  this  characteristic  carries  over  to 
the  deactivation  coefficients.  This  is  the  case  for  neat  IV’ex,  gold  foil  and  Teflon.  These 
characteristics  suggest  that  the  V-V  interactions  have  a  significant  impact  on  the  radial  (i.e., 
transverse)  distributions  of  the  populations  of  vibrationally-excited  molecules  when  the  wall 
interaction  is  strong.  In  effect,  the  V-V  processes  redistribute  the  populations  so  as  to  make 
the  observed  axial  decay  rates  slower  than  the  wall  rates.  For  weakly-interacting  surfaces,  the 
wall  effects  do  not  create  significant  radial  population  gradients  such  that  the  observed  axial 
decay  of  the  populations  track  closely  the  longitudinal  decay  of  the  populations  at  the  wall. 


45 


Table  7.  Results  of  Master  Equation  Model  fitting 


Material 

El 

kory 

V=1 

H 

■ 

11 

> 

v  =  5 

Signal 

Extrp. 

Neat 

360 

11 

kv(s‘) 

16±14% 

38±10% 

54±8% 

67±19% 

yes 

Pyrex 

^^waH) 

24 

43 

53 

77 

to 

Yva 

5.2 

13 

20 

27 

v=9 

Yv 

7.9 

15 

20 

34 

49 

-h 

340 

14 

kv(s>) 

19tll% 

38±5% 

57±5% 

75±1 1% 

no 

column 

30 

50 

70 

90 

Pyrex 

lO^Yv, 

6.4 

14.8 

25 

39 

lOVv 

11 

21 

35 

55 

99 

Gold 

351 

17 

kv(s‘) 

16.7±5% 

37±3% 

53±3% 

64±3% 

65±5% 

no 

Foil 

25 

36 

48 

58 

lOVva 

5.5 

14.5 

24.1 

34 

34 

lOVv 

8.9 

13.8 

20.9 

28 

43 

SS312 

382 

14 

kv(s-*) 

21±9% 

50±7% 

73±5% 

111±12% 

yes 

^^waB) 

47 

79 

100 

116 

to 

8.7 

28 

63 

10^ 

v=9 

loY 

25 

81 

424 

10^ 

Hi 

A13003 

375 

14 

kv(s-’) 

8.4±9% 

21±8% 

34±7% 

45±8% 

m 

^^i(wan) 

17 

33 

46 

63 

lOVva 

2.9 

7.6 

14 

20 

■ 

lOVv 

6.1 

13 

21 

34 

54 

■ 

Teflon 

334 

14 

kv(s‘) 

5.5±17% 

11±19% 

18±14% 

28±17% 

1^9 

5.7 

9.4 

14 

20 

1.7 

3.4 

6.3 

10 

■ 

BSH 

1.8 

3.0 

4.5 

6.8 

8.8 

■ 

XPS  Measurements  on  Selected  Specimens 
Examinations  of  the  surfaces  of  selected  specimens  were  carried  out  by  the  Surface 
Sciences  Laboratory  of  the  University  of  Dayton  Research  Institute/’  The  specimens 
consisted  of  witness  coupons  that  were  inserted  for  a  deactivation  experiment  in  the  catch 
tube  to  intercept  the  afterglow  flow  approximately  3  cm  after  the  laser  probe  volume  or  the 
Pyrex  specimen  tubes  were  broken  to  provide  witness  coupons  for  XPS  (x-ray  photoelectron 
spectroscopy)  examination.  Table  8  lists  the  13  specimens  that  were  examined.  The 
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objective  of  these  measurements  was  to  determine  the  degree  of  modification  of  the  surface 
constituency  by  exposure  to  active  nitrogen  either  in  the  afterglow  or  in  the  positive  column. 
This  concern  was  prompted  by  the  observation  that  the  observed  relaxation  rates  varied  over 
multiple  experiments  on  the  same  Pyrex  specimen  tube. 


Table  8.  Description  of  specimens  measured  with  XPS 


Specimen  No. 

Description  of  Specimen 

1 . tube  i 

Pyrex;  inside  near  flow  outlet  end;  Part  A 

2;  tube  i 

Pyrex;  outside  near  flow  inlet  end;  Part  B 

3;  coupon  ii 

Pyrex;  witness  in  afterglow  at  outlet  end  of  test  specimen;  Part  A 

4;  tube  ii 

Pyrex;  positive  column  exposure;  white  deposit  on  outside;  Part  B 

5;  tube  ii 

Pyrex;  positive  column  exposure;  inside  near  flow  outlet  end 

6;  tube  Hi 

Pyrex;  long  term  exposure;  inside  near  flow  outlet  end 

7;  coupon  Hi 

Pyrex;  to  check  for  metal  deposition  from  discharge  system 

8;  coupon  iv 

Pyrex;  used  with  Au  foil  specimen 

9; coupon  v 

Pyrex;  before  nitriding 

10;  coupon  vi 

Pyrex;  during  nitriding  in  afterglow;  gas  contact  (convex)  side 

11;  coupon  vi 

Pyrex;  during  nitriding  in  afterglow;  back  (concave)  side 

12;  coupon  vH 

Pyrex;  after  nitriding  in  afterglow 

13;  coupon  vHi 

Au  foU;  clipped  from  specimen  at  gas  outlet  end;  side  in  contact 
with  active  nitrogen 

Table  9  summarizes  the  results  of  the  measurements  determined  by  XPS.  This 
technique  detects  the  binding  energy  of  atoms  in  the  surface  region  to  a  depth  of  about  30 
Angstroms.  The  recorded  signal  count  rate  for  a  given  energy  is  proportional  to  the  number 
of  bonds.  The  carbon  readings  originated  from  the  material  used  to  mount  the  specimens  and-, 
therefore,  were  intrinsic  to  the  XPS  system.  There  was  no  evidence  of  carbon  on  the 
specimens.  A  good  example  of  an  interpretation  of  these  measurements  is  the  study  of  the 
inside  surface  of  specimens  Nos.  1  and  2  in  Table  8  taken  from  Pyrex  specimen  tube  i.  The 
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outlet  in  of  this  tube  had  a  1  atomic  %  of  nitrogen  bonded  to  oxygen  and  1.5  atomic  % 
bonded  to  silicon  while  the  concentration  of  the  inlet  end  was  less  than  0.3  atomic  %.  This 
result  confirmed  our  suspicions  that  the  surface  of  the  specimen  tube  was  being  modified  by 
the  afterglow.  Working  from  the  formula  and  density  of  Pyrex,  the  surface  at  the  outlet  end 
had  a  8.4%  coverage  of  nitrogen  adatoms.  A  number  of  deactivation  experiments  had  been 
run  on  this  specimen  tube  with  an  accumulated  exposure  to  the  afterglow  of  more  than  70 
hours.  These  experiments  gave  not  only  different  values  of  observed  relaxation  rates,  k^,  and 
the  corresponding  Yva^  ^^ut  different  dependencies  of  on  v.  The  results  from  the  surface 
lab  supports  the  contention  that  the  deposition  of  nitrogen  atoms  on  the  surface  affects  the 
vibrational-  energy-to-surface  exchange  mechanisms.  The  correlation  of  the  XPS  studies  with 
exposure  to  active  nitrogen  points  to  two  observations:  1)  continued  exposure  tends  to 
increase  Yu  and  2)  the  increase  of  Yva  with  v  relative  to  Yu  decreases.  As  a  general  rule,  the 
observed  deactivation  rates  k^  increased  linearly  with  v;  however,  as  the  time  of  exposure  to 
active  nitrogen  increased,  the  slope  of  k,,  with  v  decreased. 


Table  9.  Atomic  percent  of  surface  composition  of  specimens  given  in  Table  8 


No. 

Ic 

0 

Si 

N-0 

N-Si 

Na 

Cu 

Ni 

Fe 

A1 

F 

Au 

1 

■ww 

laaFa 

warn 

■sa 

2 

14.3 

60.1 

25.1 

<0.3 

<0.3 

3 

19.2 

51.0 

19.5 

— 

,0.5 

— 

5.8 

5.0 

4 

36.6 

38.5 

8.2 

— 

0.7 

15.7 

<0.3 

5 

wa 

62.9 

27.3 

1.1 

1.3 

6 

13.8 

60.5 

24.3 

1.2 

<0.3 

BH' 

21.9 

55.1 

18.9 

— 

<0.5 

1.0 

2.5 

8 

26.1 

46.2 

17.4 

— 

— 

5.8 

2.3 

— 

— 

wBSM 

— 

2.2 

9 

14.7 

58.6 

22.1 

— 

<0.6 

— 

— 

1.1 

<0.2 

Bl 

— 

0.3 

10 

11.5 

52.3 

26.1 

1.8 

— 

— 

<0.2 

<0.3 

<0.5 

— 

— 

11 

27.1 

44.9 

27.0 

— 

<0.6 

14.8 

— 

<0.2 

— 

0.8 

— 

12 

53.5 

9.5 

2.0 

1.0 

WSM 

— 

— 

— 

0.8 

■Hi! 

lUrl 

— 

BfMIM 

ERIS 

16,9 
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CONCLUSIONS 


The  results  presented  in  this  work  represent  the  first  extensive  study  of  deactivation 
of  vibrationally-excited  molecules  for  the  v  =  1,  2,  3,  4  and  5  hot  bands  on  a  number  of 
different  surfaces  correlated  with  surface  history  and  properties.  Moreover,  this  work 
includes  the  first  attempt  to  fit  a  relatively  complete  set  of  rate  equations  for  the  prevalent 
processes  to  such  data  in  order  to  extract  the  true  surface  deactivation  coefficients.  The 
values  reported  here  for  Pyrex  are  consistently  lower  for  the  v  =  1  band  by  a  factor  of  about 
two  than  the  values  published  by  Black  et  al.*  which  can,  at  least  in  part,  be  attributed  to  their 
lower  temperature  of  300  K.  This  is  supported  by  the  temperature  behavior  of  the  results 
obtained  with  the  gold  foil  wherein  the  y,,  decreases  by  almost  a  factor  of  four  for  a  change 
in  temperature  of  56  K.  The  results  obtained  for  metals  were  qualitatively  the  same  as 
reported  by  Black  et  al.  although  our  values  for  aluminum  were  somewhat  lower,  which  could 
be  due  to  different  alloy  used  here  and  our  higher  temperatures.  There  is  a  considerable 
discrepancy  between  their  value  for  Teflon  and  ours.  With  the  exception  of  titanium.  Teflon 
gave  the  lowest  yj^  values  while  Black  et  al.  found  a  value  for  Teflon  equal  to  that  for  Pyrex. 
The  results  obtained  here  for  titanium  were  quite  surprising  and  represent  a  significant  new 
result.  The  obvious  implication  is  that  a  container  of  vibrationally-excited  nitrogen  made  out 
of  titanium  could  provide  some  significant  degree  of  storage  of  this  form  of  energy. 

These  measurements  were  very  tedious  because  of  the  very  large  need  to  maintain  a 
stable  discharge-CARS  system  for  periods  extending  often  many  hours  and  because  of  the 
extensive  corrections,  normalizations  and  unfolding  of  the  data  to  obtain  the  fractional 
populations  of  the  observed  v  states  as  a  function  of  gas  residence  time  in  the  specimen  tube. 
This  tediousness  was  exacerbated  by  the  need  to  make  precise  saturation  measurements  due 
to  SRS  for  each  major  series  of  runs  that  could  be  correlated  accurately  with  the  deactivation 
measurements.  As  it  turned  out  the,  the  effect  of  the  resulting  corrections  on  the  observed" 
deactivation  rates  was  not  large,  typically  on  the  order  of  10%  or  less.  For  a  given  specimen, 
good  precision  was  usually  achieved  and  the  gold  data  showed  an  uniform  decrease  in  y^ 
with  increasing  temperature  as  would  be  expected  for  physical  adsorption. 
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A  significant  result  was  the  behavior  of  the  observed  deactivation  rates  as  a  function 
of  the  vibrational  state,  v.  It  is  typically  assumed  in  theoretical  analyses  that  the  deactivation 
coefficient,  increases  with  v  similar  to  the  increase  of  the  VT  rate  coefficient/”  For  v  <  5, 
this  rate  is  roughly  proportional  to  v.  The  observed  deactivation  rates  reported  here  do  show 
approximate  proportional  dependence  on  v  in  most  cases.  Examples  are  given  in  Fig.  7  for 
V  <  5.  Moreover,  the  deactivation  coefficients,  Yv  given  in  Table  7  (discounting  the  special 
case  of  measurements  in  the  positive  column)  that  were  derived  from  the  fitting  of  the  Master 
Equation  Model  show  only  small  deviations  from  proportionality  with  the  glaring  exception 
of  the  SS  312  results.  These  results  also  show  that  neat  Pyrex,  gold  foil  and  aluminum  all 
follow  nearly  the  same  proportional  dependence.  This  is  consistent  with  the  fact  that  two  of 
these  materials  have  very  low  reactive  oxide  surfaces  with  gold,  of  course,  being  non-reactive. 
In  these  cases,  one  expects  that  the  same  deactivation  process,  namely  physical  adsorption, 
would  be  dominant,  and  thus  show  similar  dependancies  of  the  deactivation  coefficient  on  v. 
This  result  is  all  the  more  satisfying  since  there  was  little  effort  to  control  the  surface 
roughness  of  the  specimens,  although  these  three  specimens  had  typical  smooth-appearing 
surfaces. 
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Figure  7.  Examples  of  the  linear  dependence  of  the  observed  deactivation  rates  on  v 
up  to  V  =  4  for  the  data  given  in  Figs.  4,  5  and  6 
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